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Abstract. We present new HST Space Telescope Imaging Spectrograph observations of the nearby radio galaxy NGC 5128 
(Centaurus A). The bright emission line with longest wavelength accessible from HST, [S DTJ/19533A, was used to study the 
kinematics of the ionized gas in the nuclear region with a 0'.' 1 spatial resolution. The STIS data were analized in conjunction 
with the ground-based near-infrared Very Large Telescope ISAAC spectra used by M arconi et alJ 1200 ll) to infer the presence 
of a supermassive black hole and measure its mass. The two sets of data have spatial resolutions differing by almost a factor of 
five but provide independent and consistent measures of the BH mass, which are in agreement with our previous estimate based 
on the ISAAC data alone. The gas kinematical analysis provides a mass of M BH = (1.1 ± 0.1) x 10 8 M o for an assumed disk 
inclination of i = 25 deg or Mbh = (6.5 + 0.7) x 10 7 M o for i = 35 deg, the largest ( value allowed by the data. We performed a 
detailed analysis of the effects on Mbh of the intrinsic surface brightness distribution of the emission line, a crucial ingredient in 
the gas kinematical analysis. We estimate that the associated systematic errors are no larger than 0.08 in log M BH , comparable 
with statistical errors and indicating that the method is robust. However, the intrinsic surface brightness distribution has a large 
impact on the value of the gas velocity dispersion. A mismatch between the observed and model velocity dispersion is not 
necessarily an indication of non-circular motions or kinematically hot gas, but is as easily due to an inaccurate computation 
arising from too course a model grid, or the adoption of an intrinsic brightness distribution which is too smooth. The observed 
velocity dispersion in our spectra can be matched with a circularly rotating disk and also the observed line profiles and the 
higher order moments in the Hermite expansion of the line profiles, h 3 and h 4 , are consistent with emission from such a disk. 
To our knowledge, Centaurus A is the first external galaxy for which reliable BH mass measurements from gas and stellar 
dynamics are available and, as in the case of the Galactic Center, the Mbh gas kinematical estimate is in good agreement with 
that from stellar dynamics. The BH mass in Centaurus A is in excellent agreement with the correlation with infrared luminosity 
and mass of the host spheroid but is a factor ~ 2 - 4 above the one with the stellar velocity dispersion. But this disagreement 
is not large if one takes into account the intrinsic scatter of the Mbh - o"c correlation. Finally, the high HST spatial resolution 
allows us to constrain the size of any cluster of dark objects alternative to a BH to r. < C'035 (- 0.6 pc). Thus Centaurus A 
ranks among the best cases for supermassive Black Holes in galactic nuclei. 
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1. Introduction 

One of the drivers of current astrophysical research is under- 
standing the birth and evolution of galaxies but it is becom- 
ing increasingly clear that the central supermassive black hole 
plays a major role in the evolution of the host galaxy. 

During the last few years it has been realized that most, if 
not all, nearby luminous galaxies host a supermassive black 
hole (BH) in thei r nuclei with m asses in the 10 6 - 1O 1O M 
range (e.g. lFerrarese & Fordl2005l and references therein). The 
BH mass (Mbh) is closely related with mass or luminos- 
ity of the host spheroid (L, „ i,: e.g. iKormendv & Richstonel 
1 19951 IMarconi & Hunt! 120031 Mitring & Rixll2004l) and with 
the stellar velocity d ispersion (cr: iFerrarese & Merrittl [2000 ; 
iGebhardt et alJEoOol) . These facts indicate that the formation 
of a massive BH is an essential ingredient in the process of 
galaxy formation and that there is a tight relation between host 
galaxy and the central BH. Like for t he X-ray background, 
BHs are relics of p ast AGN activity JYu & Tremainei r2002: 
IMarconi et al.ll20"04T) . and it is being realized that an accret- 
ing supermassive BH has a major impact on the star forma- 
tion rate in the host galaxy: the feedback from the accret- 
ing BH (i.e. from the AGN) is thus responsible for setting 
the c lose relations between BH and host galaxy p roperties 
(e.g. iMenci et al .1120031 loranatn et aljEool lOi Matteo et al J 
120051) . 

BHs are detected and their mass directly measured in 
nearby galaxies using gas or stars as tracers of the kinematics 
in the nuclear region. The gas kinematical method has the ad- 
vantage of being relatively simple, of requiring relatively short 
observation times and being applicable to AGNs. On the con- 
trary, non gravitational or non circular motions might bias or 
completely invalidate the method. Moreover, detectable emis- 
sion lines are not always present in galactic nuclei. The stel- 
lar dynamical method has the advantage that star motions are 
always gravitational and that stars are present in all galactic 
nuclei. However, long observation times are required to obtain 
high quality observations and stellar dynami cal models are very 
complex leading to potentia l indeterminacy JValluri et al.l2 004: 
ICretton & Emselle ml2004 . 

The weakest points of both methods - assumption of cir- 
cular gravitational motions for gas kinematics, complexity of 
modeling for stellar dynamics - have shed some doubts on 
the reliability of BH mass estimates (e.g.lTre maine et al. 2002; 
Cappellari et al.l2002tlverdoes Kleiin et alJ2002HValluri etall 
2004; Cretton & Emsellem 2004). This is a very important is- 
sue because the correlations between BH mass and host galax- 
ies properties, which are based mostly on Mbh from gas kine- 
matic and stellar dynamics, are pivotal in understanding the co- 
evolution of galaxies and BHs over cosmic time. Securing the 
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correlations means establishing without any doubt the reliabil- 
ity of gas kinematics and stellar dynamics Mbh measurements. 

Another important open issue is that the so-called super- 
massive BHs are in reality massive dark objects, i.e. clumps of 
dark mass unresolved at the spatial resolution of current ob- 
servations. Only in the case of our own galaxy, we can safely 
state that the massive dark object is a supermassive black hole. 
In all other cases the massive dark object could very well be a 
cluster of dark stellar remnants (stellar mass black holes, neu- 
tron stars, hot jupiters etc.). Such clusters, due to friction and 
collisions, are doomed to collapse to a supermassive BH but 
times cales can be longer than the age of the universe (Maoz 
119981) . In order to exclude that a massive dark object is a clus- 
ter of dark objects one needs to constrain its size such that its 
life time is much shorter than the a ge of the universe, making 
a BH a very likely alternative (e.g. [van der Marel et aflll997t 
Macchett oet al!l997llSchodel et al.l2003l) . 

Centaurus A (NGC5128) is a giant elliptical galaxy and 
hosts a powerful radio source and an AGN reveal ed by 
the p resence of a powerful radio and X-ray jet (see |lsr; 
1 for a review; see also iTingav e t al. 200 1[ | Grandi et 
lHardcastle etail l2003t lEvans et all 120041 and refer 
ences therein). The AGN in Centaurus A is l ow luminos- 
ity and is consistent with a m isaligned BL-Lac dCapetti et alJ 
2000; Chiabe rge etall l200lh . The most striking feature of 
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Centaurus A is the dust lane, rich in molecular 


ind ionized gas 


(Ouillen et al. 


1993, 1992; 


Leeuw et alJ 


2002; 


Karovska et al. 


2003; Mirabel et al. 1999), which crosses the whole galaxy hid- 



ing the nuclear region under at least 7 magnitudes of extinction 
in V JSchreier et alJI 19961: IMarconi et alJ l2000). A recent esti- 
mate of the dis tance of Centaurus A places the galaxy at ~ 3.8 
Mpc (lReikubjf 2004) but in this paper we ado pt D = 3.5 Mpc 
(1" a 17 pc. lHui et alJl993llSoria et alJl996l) to be consistent 
with our previous papers. 

ISchreier et alJ Jl998l) found evidence for a disk-like fea- 
ture in ionized gas disk from Paa HST/NICMOS images. 
IMarconi et alJ J200ll) used near infrared spectroscopy to per- 
form a gas kinematical analysis of such feature finding ev- 
idence for a BH with mass M BH = 2+]° x 10 8 M o . The 
large error associated with Mbh takes into account the de- 
genera cy due to the unknow n inclination of the gas disk for 
which IMarconi etall J200ll) could only assume i > 15deg. 
Observations were performed with a spatial resolution of about 
ft' 5 (- 8.5 pc). Even at that moderate spatial resolution the 
BH sphere of infl uence is well resolved: with cr = 138 kms 1 
JSilge et al .1120051) the radius of the BH sphere of influence is 
'"bh = GMbh/ct 2 = 1.3" . Centaurus A is thus the perfect 
benchmark to assess the validity of the gas kinematical method 
because the results will be little affected by spatial resolution or 
by uncertainties on the stellar contribution to the gravitational 
potential. The HST/STIS observations presented here have a 
spatial resolution which is a factor of 5 better than our previous 
ground-based observations allowing a unique test on the reli- 
ability of the gas kinematical method. A comparison between 
rotation curves taken at ground based and HST resolution was 
performed in Cygnus A and showed that after allowing for the 
different instrumental setups and spatial resolution, the best fit 
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ground based m odel could nicely repro duce the observed STIS 
rotation curves (Tadhunter et al. 2003). 

Recently Silg e et alJ (|2005) have presented a stellar dy- 
namical measurement of the BH mass in Centaurus A. To our 
knowledge this is the first time that a BH mass can be reli- 
ably measured in an external galaxy both with gas kinemat- 
ics and stellar dynamics and the agreement of the two mea- 
surement s can greatly strengthen the reliability of both meth- 
ods. |Caj)i>elia^euil presented gas and stellar dynam- 
ical measurements of the BH mass in IC1459 finding a dis- 
agreement of almost a factor 10 between the two. However in 
IC1459 the rotation curve of the gas is not consistent with that 
of a rotating disk and the M BH measurement cannot be consid- 
ered reliable. 

Last, but not least, the proximity of Centaurus A, combined 
with HST resolution will allow us to put a tight constraint on 
the size of any nuclear dark cluster with the possibility of estab- 
lishing whether the putative BH is really a BH and not a cluster 
of dark stellar remnants. 

In summary, Centaurus A is an ideal benchmark for direct 
BH mass measurements due to (i) its proximity, (ii) the well- 
resolved BH sphere of influence, (iii) the availability of data at 
different spatial resolution both of which resolve the BH sphere 
of influence and (iii) the availability of both gas and stellar dy- 
namical Mbh measurement. 

In Sec. |2 we present the HST/STIS observations and data 
reduction, and the re-analysis of our old VLT/ISAAC data. In 
Sec. |3 we describe the kinematical analysis, including the use 
of the Hermite expansion of the line profiles, and we present 
the kinematical parameters measured from all data. In Sec. |4] 
we describe the modeling of the kinematical data and how it is 
possible to obtain Mbh measurements from STIS and ISAAC 
data independently. In Sec. [5] we discuss the results from the 
modeling of STIS and ISAAC data, we evaluate the influence 
of the intrinsic surface brightness distribution of the disk on the 
derived Mbh value, and we present the problems which might 
be encountered when comparing observed and model velocity 
dispersions of the gas, the commonly adopted test of the re- 
liability of the gas kinematical analysis. We compare our Mbh 
value with the stellar dynamical measurement, we discuss Mbh 
within the context of the correlation with host galaxy parame- 
ters and we then estimate an upper limit on the size of any dark 
cluster alternative to a BH. Finally, in Sec.|6]we summarize our 
work and present our conclusions. 

2. Observations and Data Reduction 

2.1. STIS observations 

Centaurus A was observed with STIS on the HST in 2000 
May 9 and July 1. The nucleus location in Centaurus A is 
not well defined in optical (< 1 /im) images thus telescope 
acquisition was performed on a foreground star ~ 10" away 
from the galaxy nucleus in order to register the coordinate 
system of the telescope (Aa, Ad = -9.3", -3.3" in Fig. [^). 
During each visit, a ~ 5" x 5" acquisition image was ob- 
tained with the F28X50LP filter on the foreground star which 
was subsequently centered and re-imaged following the ACQ 



Table 1. Log of observations. 



ID 


Date 


Instr. 


PA 


'R J 


Res b 








[deg] 




["] 


VIS1 


May 9, 2000 


STIS 


-5.0 


1000 


0.1 


VIS2 


July 1, 2000 


STIS 


43.0 


1000 


0.1 


PA1 


July 21, 1999 


ISAAC 


32.5 


10000 


0.5 


PA2 


July 22, 1999 


ISAAC 


-44.5 


10000 


0.5 


PA3 


July 22, 1999 


ISAAC 


83.5 


10000 


0.5 



a Spectral resolution. 



b Spatial resolution: FWHM of Point Spread Function (PSF) or seeing. 

and ACQ/PEAK procedures. The exposure times of the acqui- 
sition images were 3 and 10 s respectively. After registering 
the coordinates system on the star, the telescope was subse- 
quently pointed on the expected position of the Centaurus A 
nu cleus derived from the analysis of WFPC2 images performed 
bv lMarconi etall l|2000). The observational strategy then con- 
sisted in obtaining an image in the F28X50LP filter to check the 
telescope pointing and slit locations with respect to the nucleus. 
Spectra were then obtained with the 52X0. 1 aperture (07 1 slit 
width) centered on the expected nucleus position. Given the 
different roll angle of the telescope during the two visits, the 
slits were located at PA -5 deg (hereafter VIS1) and 43 deg 
(hereafter VIS2). The slit positions are overlaid on the STIS 
image in Fig. At each slit position we obtained three spec- 
tra with the G750L grating centered at 7751 A, in order to ob- 
serve the [S m]M 907 1,9533 A doublet. The second and third 
were spectra shifted along the slit by 075 and 1" in order to 
remove cosmic-ray hits and hot pixels. The spectra were ob- 
tained with the 07 1 slit and no binning of the detector pixels, 
yielding a spatial scale of 070507/pix along the slit, a disper- 
sion per pixel of AA = 4.882 A and a spectral resolution of 
% = A/(2AA) = 1000 at 9500A. The log of the observations 
is shown in Tabled Th e raw STIS images were processed 
with the calstis pipeline dKim Ouiiano et ai]l2003l) using the 
best reference files available as of December 2003. The raw 
STIS spectra were processed with standard pipeline tasks but 
following a different procedure to correct the fringing exhib- 
ited by the STIS CCD in the far red (> 7500A). In order to cor- 
rect for fringing, we obtained contemporaneous flat fields with 
the same setup as the scienti fic observations. We then followed 
the prescriptions of Goudfrooii & Christensenl (11998). The on- 
orbit flats were normalized using the task normspf lat in the 
stsdas.hst_calib.stis package of IRAF 1 . Spectra were 
first prepared with prepspec and defringed with defringe. 
We did not use the task mkfringeflat to rescale and shift 
the flat fringes to those observed in the object because the 
signal-to-noise ratio on the continuum was very poor. However, 
using two flats obtained in different spacecraft orbits during 
each visit, we verified that variations of the fringe positions, 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. Left: slit positions overlaid on the Centaurus A STIS image. VIS1 and VIS2 identify the STIS slit positions and PA1, PA2 
and PA3 the ISAAC ones. The 0,0 position is the target position assumed by STIS and derived by the location of the bright star 
on the right (located at -9", -3"). Right: Zoom on the central 1" x 1". 
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Fig. 2. Average STIS spectrum from VIS 1 and VIS2 extracted from aff.'l x 076 aperture centered on the assumed nucleus 
position. 



if any, do not affect the final rotation curves. The defringed 
spectra were then wavelength calibrated (wavecal) and cor- 
rected for 2D distorsion (x2d). The spectra taken at different 
slit positions were then raligned using as reference the [S III] 
line and combined rejecting deviant pixels. The expected ac- 
curacy of the wavelength calibration is 0.1 - 0.3 pix within 



a single exposure and 0.2 - 0.5 pix among different exposures 
(IKim Ouiiano et al. 2003) which converts into ~ 15-45 km s 
(relative) and ~ 30 -70 km s" 1 (absolute). The relative error on 
the wavelength calibration is negligible for the data presented 
here because our analysis is restricted to the small detector re- 
gion including [S III] (AA < 300A). From an analysis of the arc 
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STISVIS1 STISVIS2 




Fig. 3. Examples of the fits of the line profiles for the STIS observations of VIS1 and VIS2. The dashed and dot-dashed lines 
identify the two components of the fit, called 'main' and 'blue'. For every slit position the blue component has the same velocity 
and width at each row and their values are determined in the fit of the overall nuclear spectrum shown in top panels. The numbers 
in the upper right corners of each panel are the row and range of rows where the fit was performed. 



exposures it is found that the lamp lines at ~ 9000 - 9500A 
have a FWHM of ~ 2.2 pix consistently w ith the expected 
spectral resolution JKim Ouiiano et af]l2003h . Therefore, the 
degradation of the spectral resolution at ~ 9000 - 9500A pro- 
duced by the grating corresponds to a velocity dispersion of 
cr = 140 km s -1 . When also the degradation introduced by the 
slit width is considered, the final instrumental resolution cor- 
responds to <x = 160kms~', i.e. the width of the arc lines at 
~ 9000 - 9500A. 

The target acquisition strategy followed in these observa- 
tions ensures that the STIS slits are centered on the same point, 
regardless of its distance from the Centaurus A nucleus. This 
can be verified from the images obtained just after the telescope 
offset from the acquisition star and before the spectra. The aper- 



ture center coordinates on the images are given by the HOSIAX 
and HOSIAY keywords. We used these nominal aperture centers 
as pivots for rotation and realignment of the images taken in 
the 2 different visits. By using 9 point sources common to both 
images we then verified that the images are aligned to better 
than 0.1 and 0.25 pixels rms in RA and DEC, i.e. 07006 and 
07014 respectively. This is then the accuracy with which the 
STIS slits are centered on the same point which is shown as a 
filled square in Fig.^). 

2.2. ISAAC Observations 

We complement the new HST STIS Observations with the J- 
band, ground-based spectroscopic data already presented and 
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Fig. 4. From top to bottom: Average emission line surface brightnesses (units of 10 16 ergs 1 cm 2 arcsec 2 ), velocities, and 
velocity dispersion observed along the slit in STIS observations. Columns 1 to and 2 refers to VIS1 and VIS2, respectively. 
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Fig. 5. Examples of the fits of the line profiles for the ISAAC observations of PA1, PA2 and PA3. The solid lines represents the 
case of the Hermite expansion while the dashed line is the simple gaussian fit. The line at shorter wavelengths is [Fe II] while at 
longer wavelengths the strongest line is Pa/?. All emission lines have the same average velocity, velocity dispersion and Hermite 
coefficients 

. At the blue side of [Fe II] a faint line is visible, this is [S IX]. The line at the blue side of Pa/? is another [Fe II]. The dot dash short line mark 

the expected position of the blue component found in the [S III] line. 



analyzed bv iMarconi et alJ {2001). Briefly, the data were ob- 
tained in 1999 using ISAAC at the ESO VLT and consist of 
medium resolution spectra obtained with the 0/3 slit and the 
grating centered at A = 1 .27 fim. The dispersion is 0.6 A/pix, 



yelding a resolving power of 10000. Seeing during the obser- 
vations was in the range 0'.'4-0'.'6, with photometric conditions. 
The observation procedure was to obtain an acquisition im- 
age in the K band and then center the slit on the prominent K 
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Fig. 6. From top to bottom: Average emission line fluxes, velocities, velocity dispersion, hi and h\ observed along the slit in 
ISAAC observations. Columns 1 to 3 refers to PA1, PA2 and PA3, respectively. The kinematical quantities were derived from a 
simultaneous fit of Pa/? and [Fe II]. 
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band peak in the nuclear region of Centaurus A. Spectra were 
then obtained with the slits placed at three different position an- 
gles, summarized in Tabled More details on observ ations and 
data reduction can be found in M arconi et alJ J200ll) . We have 
re-re duced the ISAA C spectra following the same strategy as 
iMarconi et alJbOOll) but with a difference: before coadding all 
the wavelength calibrated frames, we have traced the position 
of the nuclear continuum peak along the dispersion direction 
and used this trace to rectify the dispersion direction. This en- 
sures that there is no variation of the location of the continuum 
peak along the slit in the wavelength range covered by the ob- 
servations. 

3. Results 

3.1. STIS Line Kinematics 

Fig.|2shows the average STIS spectrum from VIS1 and VIS2 
extracted from a CC.'l x (X'6 aperture centered on the assumed 
nucleus position. The lines Ha, [Nil], [SII], [Oil], [S III] 
and [CI] are clearly detected and their relative strenght are 
an indication of the strong reddening which affec ts the nu- 
clear region of Centaurus A (A v > 7 mag e.g. JSchreier et alJ 
Il996l IMarconi et alJEoOol) : typically, Ha+[NII] is stronger 
than [S III] while here [S III] /Ha + [N II] > 10. Therefore, only 
[S III]/t9500A has enough signal-to-noise ratio (hereafter SNR) 
for kinematical analysis. We thus focus only on the 9300-9800 
A spectral region. 

After selecting the spectral region, we subtracted the con- 
tinuum by fitting a linear polynomial row by row along the dis- 
persion direction. The continuum subtracted lines were then fit- 
ted row by row with gaussian functions using the task specfit in 
the IRAF stsdas.contrib package. When the signal-to-noise ra- 
tio (SNR) was insufficient (faint line) the fitting was improved 
by co-adding two or more pixels along the slit direction. 

Fitting single gaussians is acceptable for most of the po- 
sitions along the slit but it does not produce good results for 
the points in the nuclear region where a careful analysis of 
the line profiles shows that they are persistently asymmetric 
with the presence of a blue wing (Fig. [3}. A fit row by row 
with two gaussian components, the main component and the 
blue one, shows that, within the large uncertainties, the "blue 
wing" has always the same velocity and width. We have thus 
deblended the "blue" component in the spectrum obtained by 
co-adding the central 7 rows. The velocity and width of the blue 
component (v = — 240kms _I and <x = 360kms~ I for VIS1, 
v = -320kms _1 and <x = 430kms~ 1 for VIS2) were then held 
fixed in the row by row fit significantly improving the fitting of 
the line profiles (Fig. QJ. A similar blue wing was detected and 
deblended in Ha+[NII] in NGC 4041 (IMarconi etalj E003a) 
and, as in that paper, the present data does not allow to identify 
the nature of the blue wing. We remark that the final rotation 
curve does not depend on the exact values of v and cr used to 
parameterize the blue wing since that component is faint com- 
pared to the main one. However, the rotation curves obtained 
without taking into account the blue wing are not consistent 
with a circularly ro tating disk, like in the case of NGC4041 
iMarconi et al.l2003al). 



The fitting procedure provides line-of-sight velocities, ve- 
locity dispersions and surface brightnesses along each slit for 
the [S III]/19500A emission line which are shown in Fig.|4] 

The observed line surface brightness in VIS 1 and VIS2 are 
strongly peaked at the expec ted nucleus position s, similarly 
to what was observed in Paa JSchreieretalJll998l) and [Fell] 
jMarconi et alJl2000l l200lh . The SNR of the data does not al- 
low the detection of extended emission beyond 1" from the nu- 
cleus. The rotation curve in VIS1 is "S-shaped" as expected in 
the case of a circularly rotating disk observed with a slit close 
to its major axis. In VIS2 the S-shaped pattern seems superim- 
posed on a constant velocity gradient. The amplitude of the ro- 
tation curves is ~ 300 km s" 1 (VIS1) and ~ 200 km s 1 (VIS2) 
suggesting that the VIS 1 slit is closer to the orientation of the 
disk line of nodes. The behaviour of the velocity dispersion is 
similar to what is qualitatively expected from a rotating disk, 
i.e. a constant <x which increases toward the nucleus location 
reaching a peak value of ~ 350 km s _1 . 

3.2. ISAAC Line Kinematics 

The fit with one or two gaussians - in case of a significant 
blue wing - produces good results for the STIS data because 
of the moderately low SNR and spectral resolution (~ 1000). 
However, the ISAAC data have a much higher SNR and spec- 
tral resolution and the line profiles present a significant struc- 
ture in the nuclear re gion similar to the exp ected one in case of 
unresolved rotation. IMarconi et alJ J200ll) estimated kinemat- 
ical parameters by fitting single gaussians to [Fell] and Pa/3 
lines. While the fitting was acceptable for the purposes of that 
paper, it had significant residuals. In order to estimate accu- 
rately the average velocity and velocity dispersion, i.e. the same 
parameters computed by our model (see below), the residuals 
should not have any structure. When a simple gaussian is a too 
simple model for the observed line profile, a possibility is to fit 
2 gaussians, but the fitting is unstable and gaussian parameters 
are strongly correlated. Therefore we decided to model the line 
pr ofile with the expansion in He rmite polynomial introduced 
bv Ivan der Marel & Franxl ill 9931) coupled with the penalized 
pixel fitting technique JMerrittlll997t Cappella ri & EmsellerrJ 
2004). Such parameterization of the line profile coupled with 
the penalized fitting technique has been used to describe the 
Line of Sight Velocity distribution of stellar kinematical data 
but it is easily applicable to the analysis of gas kinematical data. 
Of course, this refinement is necessary only because of the sig- 
nal to noise and spectral resolution of the data is moderately 
high, otherwise a single gaussian woul d be more than adeq uate 
to des cribe the data. Briefly, following Ivan der Marel & Franxl 
we parameterize the line profile L(v) as: 

L(v) = ya(w) [1 + h 3 H 3 (w) + h 4 H 4 (w)] (1) 

where 



a(w) = ^-e-^- w2 (3) 
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Fig. 7. Comparison of averag e velocities from [Fe I II and Pa/3 
emission lines with those from lMarconi et aDfcOOll) . 



and Ht,(w), H^iw) are the third and fourth order 
Her mite polynomials with t he normalization adopted by 
Ivan der Marel & Franxl dl993h . vo and <xo are the central 
velocity and velocity dispersion of the gaussian component 
and correspond to the average velocity and velocity dispersion 
of the line profile when it is a perfect gaussian. and /14 
quantify the deviations from the gaussian line profile and 
/13 = /14 = corresponds to the perfect gaussian case. Though 
the Hermite expansion can be extended to orders higher than 
4 th , we did not find significant improvements in the profile 
fitting, given the SNR of the data. The above line profile was 
fitted by minimizing the penalized x 1 where 



4 



-■X 2 [l+P 2 (hl + h 2 4 )] 



(4) 



see [baprjellari & Emsel leml J2004h for more details, x 2 is the 
classical Chi-squared and f3 is a "bias" factor which favours 
best fits with low h.3 and h\ values. We found that ft — 1 is 
adequate to the spectral resolution and SNR of our data. We 




400- 
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Position along the slit [arcsec] 

Fig. 8. Comparison of average velocities from STIS VIS2 and 
ISAAC PA 1. 



investigated comprehensively the possibility that lines did not 
share common kinematics but we did not find evidence for this. 
On this basis we proceded the subsequent analysis assuming 
that all lines in the spectrum have the same average velocity, 
velocity dispersion, hi and h\. 

Examples of the fits of the line profiles for the ISAAC ob- 
servations of PA1, PA2 and PA3 are shown in Fig [5] where we 
also compare the Hermite expansion with the case of the pure 
gaussian fit (dashed line). 

The fitting procedure provides line-of-sight velocities, ve- 
locity dispersions and surface brightnesses along each slit for 
the [Fe II] and Pa/? emission lines which are shown in Fig.|6| In 
Fig.0we comp are our new determina tion of the average veloc- 
ity with that by M arconi et alJ J200ll) . The new rotation curves 
does not differ significantly from the old ones apart from being 
less noisy. I SAAC kinematical da ta are described and discussed 
at length in iMarconi et.al l fcOOlh and we refer to that paper for 
more details. Here we concentrate exclusively in the nuclear 
region, i.e. < 2" from the position of the K band peak. In the 
nuclear region the ISAAC data exhibit the same behaviour as 
the STIS ones, i.e. strongly peaked surface brightness profiles, 
S-shaped rotation curves with amplitudes that vary with slit PA 
(~ ISOkms- 1 at PA1, ~ 200 km s" 1 at PA2, and ~ SOkms" 1 
at PA3). The observed velocity dispersion is characterized by a 
well defined constant plateau, expecially in PA1, and a smoth 
increase toward the nucleus location where it reaches an am- 
plitude of ~ 250kms _1 at all slit orientations, hi and Ka are 
significantly different from only in within 1" of the nucleus 
where the SNR of the data is higher. In all other locations they 
have been fixed to 0. 

In Fig.|8]we compare the rotation curve from STIS at VIS2 
and the rotation curves from ISAAC at PA1. This comparison 
is not affected by geometrical effects because the two slits have 
a similar orientation, differing only by 5 deg. The ISAAC data 
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Fig. 9. Top, from left to right: observed NICMOS F222M image, dereddened F222M image, derived E(B-V) map. Bottom, from 
left to right: observed 2MASS K band image, dereddened K image, derived E(B-V) map. 



are characterized by a much higher SNR and a by a smoother 
rotation curve. The effect of the different spatial resolution and 
slit width is clearly that of smoothing away the high velocity 
increase observed in the S TIS data. This effect is ex pected and 
was discussed at length bv lMacchetto et alJ dl997l) . The STIS 
rotation curve is thus in good agreement with the ISAAC data 
apart for the two couples of points +077 away from the center. 
It is difficult to establish whether the discrepancy is due to low 
SNR of the data and/or to the presence of non circular motions 
in the [S III] emitting medium. Since we have verified that the 
inclusion of those points does not change the final results apart 
for incresing x 2 , we have decided not to consider them in the 
following analysis. 

4. Modeling the rotation curves 

The aim of our new HST/STIS observations is to confirm at 
higher angular resoluti on the detection of t he supermassive 
BH in Centaurus A bv lMarconi etalJ l l200ll) . Both STIS and 
ISAAC data are analized to verify if they provide consistent 
results. 



To model the kinemati cal data we follow the procedure 
first described by Mac chetto et al.ldl997l) and subsequently re- 
fined by several authors (Ivan der Marel & van den BoschTl 998 : 
iMarconi etai1l200lt fearth et al]l200lt IMarconi et al.ll20Q3al). 
Our m odeling code, described in detail bv HMarcorn^tall 
(2003a|), was used to fit the observed rotation curves. Very 
briefly the code computes the rotation curves of the gas assum- 
ing that the gas is rotating in circular orbits within a thin disk in 
principal plane of the galaxy potential. We neglect any hydro- 
dynamical effect like gas pressure. The gravitational potential 
is made of two components: the stellar potential, characterized 
by its mass-to-light ratio M/L, and a dark mass concentration 
(the putative black hole), spatially unresolved at HST+STIS 
resolution and characterized by its total mass Mbh- In comput- 
ing the rotation curves we take into account the finite spatial 
resolution of the observations, the intrinsic surface brightness 
distribution of the emission lines (hereafter ISBD) and we in- 
tegrate over the slit and pixel area. The free parameters charac- 
terizing the best fitting model are found by standard x 2 mini- 
mization. 
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Fig. 10. Results of ellipse fitting on dereddened NICMOS F222M (left) and 2MASS K images (right). The top panels show the 
observed surface brightness as a function of radius. Middle panels and lower panels show ellipticity and position angle of the 
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In Sec. 14. II we determine the radial light profile of the stel- 
lar component which, multiplied, by M/L directly provide the 
enclosed mass at a given distance from the nucleus. In Sec. 14.21 
we determine the ISBD which is the weight for averaging the 
kinematical quantities. Finally, in Sec. 14.31 we present the re- 
sults of the kinematical fitting obtained using the stellar mass 
distribution and intrinsic surface brightness previously deter- 
mined. In the same section we also present a new way to com- 
pute the kinematical model quantities in order to avoid subsam- 
pling problems. 



potential axial ratio q, and the inclination of its principal plane 
with respect to the line of sight, i. As these two quantities are 
related by the observed isophote ellipticity, we are left with the 
freedom of assuming differ ent galaxy inclinations to the line 
of sight. Following lvan der Marel & van den Boschl Jl998). we 
assumed an oblate spheroid density distribution parameterized 



as: 



(m 
— 
n 



1 + 1™ 

rb 



(5) 



4.1. The stellar mass distribution 

aspect 

The inversion procedure to derive the stars distribution 
from the surface brightness is not unique if the gravitational 
potential does not have a spherical symmetry, revealed by cir- 
cular isophotes. Assuming that the gravitational potential is an 
oblate spheroid, the inversion depends on the knowledge of the 



m is given by m 2 = x 2 + y 2 + z 2 /q 2 where xyz is a reference 
system with the xy plane corresponding to the principal plane 
of the potential and q is the intrinsic axial ratio. This three- 
dimensional density model is converted to an observed sur- 
face brightness distribution in the plane of the sky by integrat- 
ing along the line of sight, convolving with the Point Spread 
Function (PSF) of the telescope+instrument system and aver- 
aging over the detector pixel size. Then, the derived model light 
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Fig. 12. Upper panel: brightness profile of Centaurus A in 
the K band resulting from merging the NICMOS F222M 
and 2MASS K profiles. Surface brightness is in units of 
ergs 1 cm -2 A -1 arcsec' 2 . The solid line is the fit to the profile 
obtained from a spherical stellar density distribution with an 
added nuclear point source. The dotted and dashed lines show 
the relative contributions of point source and stellar density, re- 
spectively. 



profile can be directly compared with observed ones. A detailed 
description of the relevant formul as and of the inver sion and fit 
procedure is presented in the lMarconi et al.l {2003a). 

We reconstructed the galaxy light profile combining 
a NICMOS F222M (K band) im age obtained with the 
NIC2 camera (pixel size 07075; see ISchreier et alJll998l and 
iMarconi et"ai1l2000l for details) with a 2M ASS K band imag e 
from the Large Galaxy Atlas (pizel size 1"; Jarr ett et alJ2 003). 
The first image, with smaller pixel size and better spatial res- 
olution was used for the central regions (r < 6") while for the 
more extended emission we took advantage of the larger field 
of view of the 2MASS image. 

The nuclear region of Centaurus A is strongly reddened 
due to the well known dust lane. This effect can be seen also 
in the K band and significantly affects the profile determina- 
tion. Therefore we ap plied a redden ing correction to the im- 
ages, following M arconi et a l. (2000) where the reader should 
refer to for more details. Briefly, using available NIC2 F160W 
and 2MASS J band images of Centaurus A, we computed 
the H - K and J - K colors. Assuming as intrinsic colors 
H — K — 0.2 and J - K — 0.92 we can then estimate red- 
dening using the extinction curve bv lCardelli etalH l989) with 
Ay/E(B - V) = 3.1. Fig.[9]shows the effects of this reddening 
correction. The isophotes from the NIC2 F222M image prior to 
correction are elongated in a direction perpendicular to the dust 
lane, while after correction their circular simmetry is recovered. 
Similarly, in the 2MASS K band image prior to correction one 
can clearly see the distorsion caused by the lower edge of the 
dust lane running approximately at PA -45 deg. After correction 
the circular simmetry is well recovered in the central region. 
There is still a residual distortion due to an over correction in 
the region where the dust screen approximation fails. However 
this distortion can be easily taken into account in the ellipse 
fitting described below. 

We used the IRAF/STSDAS program ellipse to fit ellip- 
tical isophotes to the galaxy after reddening correction and the 
results are shown in Fig. 1101 

The surface brightness profiles derived from NICMOS data 
show a central unresolved source clearly identifiable by its first 
Airy ring at ~ 073. The ellipticity is small being < 0.1 if one 
excludes the inner 074 where emission is dominated by the 
central unresolved source. The position angle shows significant 
variations but it is probably affected by residual reddening. The 
surface brightness profile, ellipticity and PA show a smoother 
behaviour in the 2MASS data. The ellipticity is still small and 
increases to 0. 1 only at scales larger than those probed by our 
spectra (> 20")- 

As shown in Fig.^^the 2 dereddened profiles are in excel- 
lent agreement in the overlap region taking into account that no 
rescaling of the data was performed. 

In order to determine the stellar contribution to the mass 
density we can then safely assume a spherical simmetry since 
the isophotes, at least in the nuclear region, are close to circular. 
Moreover, we can use an observed light profile obtained by 
merging the NICMOS (r < 6") and the 2MASS ones (r > 
6"). Since the spatial resolution of the observations affects only 
the inner parts of the light profiles, the PSF to be used in the 



A. Marconi et al.: The Black Hole in Centaurus A 



13 




Fig. 14. Reference frame xy in the plane of the sky. x and y are aligned along the North and East directions, respectively. Left: 
xy characterization of the location of an elliptical isophote of the z'-th component of the ISBD. Right: xy characterization of the 
location of the j-ih slit. 
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Fig. 13. Circular velocity derived from the model fit in Fig. 1121 
and computed assuming M/L = 1 (solid line). The dashed line 
represents the rotation curve expected in the case of a BH with 
mass Mbh=1O 8 M . The small inset shows the behaviour of the 
circular velocity in the nuclear region. The radius of the BH 
sphere of influence, i.e. the radius at which the contribution of 
stars and BH to the circular velocity are equal is ^ 17 1 . 



fittin g is the one of N IC2 F222M, which we take from Tiny Tim 
jKrist&Hookll999l, . 

The best fit of the observed stellar light profile is shown in 
Fig. 1121 and is characterized by po = (1000 + 3)M Q pc~ 3 (for 
M/L = I), a = 0.81+0.01,/? = 0.74+0.01 andr,, = 3762+0701 
(~ 62 pc at the distance of Centaurus A). The lcr statistical 
errors were estimated with 100 montecarlo evaluations of the 



fit. The presence of an unresolved nuclear point source, and 
its associated Airy ring, are clearly visible in Fig. EH Indeed 
a point source with flux (3.39 + 0.01) x 10~ 15 erg s -1 cirT 2 A _1 
(K — 10.1 mag) had to be added to the extended luminosity 
distribution to provide a good fit to the brightness profile. The 
fit is excellent everywhere except for large residuals in the re- 
gion where the central unresolved emission dominates. This is 
not a serious issue since it is related to the detailed shape of 
the PSF, which is likely distorted by the reddening correction 
procedure described in Sec. 14.11 an d does not aff e ct the stel- 
lar emission. It has been shown bv lOuillen et alJ j^QOlD that 
unresolved infrared sources are found in the great majority of 
HST NICMOS images of Seyfert galaxies and that their lu- 
minosities strongly correlate with both the hard X-ray and the 
[OIII] line luminosity. This result suggests a dominant AGN 
contribution to the IR emission. Moreover, we found that the 
central unresolved s ource is significantly polarized P — 11.1% 
dCaoetti et alJ2 000) and this definitely excludes that it could be 
emission from an unresolved star cluster. 

Finally, in Fig. ^] we plot the circular velocity due to the 
stars as a function of the distance from the nucleus, com- 
puted assuming M/L = 1 (solid line). The dashed line rep- 
resents the circular velocity expected from a BH with mass 
M B h=1O 8 M . Below 1", the BH completely dominates the 
gravita tional potential, a fact already noted by iMarconi et alJ 
J200ll) . At r = 1.1", the contribution of BH and stars to the 
gravitational potential are equal. This size can be compared 
with the standard definition of the radius of the BH sphere of 
influence, tbh = GMsu/crl, where G is the gravitational con- 
stant and cr e is the effective velocity dispersion of th e stars. In 
the ca se of Centaurus A, cr e = 138 + lOkms -1 dSilge et alJ 
2005, see Sec|53J which combined with M BH =1O 8 M yelds 
tbh = 1 .3 + 0.3" in excellent agreement with the estimate from 
our analysis of the stellar potential. 



14 



A. Marconi et al.: The Black Hole in Centaurus A 




Fig. 16. Fit of the ISAAC line surface brightness distribution along the slit for 2 sample cases. 



4.2. The intrinsic line surface brightness distribution 

The observed kinematical quantities are averages over aper- 
tures defined by the slit width and the detector pixel size along 
the slit. The ISBD is thus a fundamental ingredient in the mod- 
eling of the kinematical quantities because it is the weight of 
the ave raging process (se e, e.g., the discussio n in lMarconi et alJ 
l2003al and Appendix Bl. lBarth et ail feOOll) derived the ISBD 
by deconvolving a continuum subtracted emission line image 



and were successful in reproducing the microstructure of the 
rotation curves thus reducing the overall x 2 of their best fit- 
ting model. However, apart for the microstructure of the rota- 
tion curves which has little impact on the final Mbh estimate, 
iMarconi et alJ d2003all showed that it is crucial to adopt a good 
ISBD in the central region where the velocity gradients are the 
largest. This has a strong impact on the quality of the fit but 
also on the final Mbh estimate. When the observed line surface 
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Fig. 17. Left: isophotal contour levels of the intrinsinc surface brightness distributions G2 and E4 used in Fig.[2]to reproduce 
the STIS data. The ISBDs have not been convolved with the instrumental response. Contour levels starts from -2 and go to -0.2 in 
log of the peak of the ISBD. The filled circle indicate the location of the slit centers. Right: as in the left panel but for the ISBDs 
G7 and G7b used in Fig.^]to reproduce the ISAAC data. 



Table 2. Intrinsic line surface brightness distributions for STIS 



id 


function 


( 


I 0i 


rot 


*0i 


yoi 


e, 


9i 


G2 


gauss 


1 


3685 


0.145 


0.028 


0.004 


32.5* 


0.5* 




gauss 


2 


1950 


0.231 


-0.178 


0.196 


0.00* 


1.0* 




gauss 


3 


295 


0.278 


-0.481 


-0.311 


0.00* 


1.0* 




const 


4 


90 












E4 


expo 


1 


2438 


0.066 


0.016 


0.036 


77.5* 


0.5* 




expo 


2 


3333 


0.056 


0.023 


0.006 


-12.5* 


0.5* 




expo 


3 


1948 


0.058 


-0.217 


0.103 


0.00* 


1.0* 




expo 


4 


440 


0.097 


-0.458 


-0.352 


0.00* 


1.0* 




const 


5 


83 













* fixed parameter. 



Table 3. Intrinsic line surface brightness distributions for ISAAC 



id 


function 


hi 


roi 


XQi 


yoi 


0i 


1i 


G7 


gauss 


1 40500 


0.119 


0.00* 


0.00* 


122.5* 


0.5* 




gauss 


I 591 


1.262 


0.00* 


0.00* 


32.5* 


0.5* 




gauss 


68 


3.407 


0.00* 


0.00* 


0.00* 


1.0* 


G7b 


gauss 


1 41184 


0.121 


0.106 


0.009 


12.5* 


0.5* 




gauss 


1 627 


1.190 


0.00* 


0.00* 


32.5* 


0.5* 




gauss 


64 


9.478 


0.00* 


0.00* 


0.00* 


1.0* 



* fixed parameter. 



brightness distribution along the slit is strongly peaked in the 
nuclear region and has a profile which is little different from 
that of an unresolved source, as in the present case, one cannot 
simply use a deconvolved emission line image but should try 
to reproduce the observed line surface brightness distribution 
with a parameterized intrinsic one and taking into account the 



intrumental effects. Therefore we decided to parameterize the 
ISBD with a combination of analytic functions. We will show 
in the following analysis that the final result, i.e. the BH mass, 
is very little dependent on the assumed ISBD. 

We adopt the reference system xy in the plane of the sky 
where the x and y axes are aligned with the North and East 
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directions, respectively (see Fig. I14> . The total emission line 
surface brightness at a point x, y is represented in our model by 



(6) 



where /o, is the amplitude of the /-th component which is de- 
scribed by the analytic function fi(ri/roi). ro, is a scale radius 
and r, is the distance from a simmetry center defined by 

x\ — X — XQi 

y'i = y-yoi 

x" = x'j cos 9, + y] sin ft 
y" = (-X- cos 6j + y'j sin 

il/2 



r» - [(*f) 2 + 0>f) 2 ] 



(7) 



Thus, the isophotes of the ;-th surface brightness component 
are ellipses centered in (;to;, yod, with axial ratio q/ (minor over 
major axis) and major axis aligned along a direction with po- 
sition angle PA = 9, (Fig. 1141. The adopted analytic functions 
are exponentials, gaussians and constants defined as: 

Mn) = exp(-r,/r 0/ ) 
1 



fi(n) = exp 
fi(n) = 1 {n <= r m ); (n > r m ) 



(8) 



The ISBD of Eq.gjis then convolved with the PSF of the 
system and integrated over apertures defined by the slit and 
the detector pixel size. The free parameters of the fit are de- 
termined with a x 1 minimization by comparing observed and 
model data. The position of the y'-th in the xy reference system 
is characterized by X(sy,y(jS)j, the slit center position and %y 
its position angle on the plane of the sky. We now derive the 
ISBDs for STIS and ISAAC data. 

4.2.1. STIS data 

The acquisition procedure followed during HST observations 
(i.e. centering on a bright star and slew to the expected nucleus 
location) ensures that the STIS slits are centered on the same 
point within an accuracy of a fraction of a pixel. Inspection 
of Fig. 0]then suggests that the ISBD should be described at 
least by 4 components: one at the slit center (the nuclear com- 
ponent), one offset by -072 along the slit at VIS1, one at -076 
along the slit at VIS2 and a constant base. The flattening and 
orientation of the offnuclear components can be left as free pa- 
rameters of the fit. However, in the case of a nuclear compo- 
nent with a size smaller than the spatial resolution it is safer 
to reduce the number of free parameters and therefore we have 
considered the following cases for the nuclear component: 

(1) nuclear component circularly symmetric; 

(2) nuclear component with PA = 32.5 deg and q = 0.5; 

(3) nuclear component with PA = 122.5 deg and q = 0.5; 

(4) 2 nuclear components with PA = 77.5, -12.5 deg and q = 
0.5; 

(5) 2 nuclear components with PA - 77.5, -12.5 deg and q = 
0.17; 



(6) 2 nuclear components circularly symmetric; 

(7) 2 nuclear components, one circularly symmetric and the 
other with PA - 32.5 deg and q - 0.5; 

(8) 2 nuclear components, one circularly symmetric and the 
other with PA = 122.5 deg and q = 0.5; 

(9) 2 nuclear components with PA = 32.5, 122.5 deg and q = 
0.5; 

PA - 32.5 deg is the PA of the Centaurus A jet and is the 
direction along which Pact and [Fell] emissions are elon- 
gated and q = 0.5 is their axial ratio dSchreier et al.l Il998l 
iMarconi etakl koOO). PA - 122.5 deg is the direction perpen- 
dicular to the previous one. In the case of 2 components with 
PA - 77.5,-12.5 deg we mimick an ionization cone with axis 
at PA — 32.5 deg and opening angle of 45 deg. 

For each of the above cases we have used exponential (E), 
gaussian (G) and constant (C) functions. An ISBD is then la- 
beled as, e.g, E5 which means exponential functions in case 5 
above. 

Sample results of the fitting procedure are shown in Fig. 1151 
for the ISBDs G2 and E4 and the best fit parameters for these 
ISBDs are shown in Tab.|2] Inspection of Fig.[T5lindicates that 
ISBDs G2 and E4, though different, result in very similar ob- 
served surface brightness distributions. 

4.2.2. ISAAC data 

The observational procedure followed during IS AAC observa- 
tions (i.e. centering the slit on the K band peak. liMarconi et all 
12001 1) ensures that the three slits are centered on the same point 
within a fraction of the ISAAC pixel. However it is not clear if 
the emission line peak obs erved in Pa/? and [Fe I II is coincident 
with the continuum peak. Schrei er" st al.l {l998) analyze a Paa 
image and find that the continuum subtracted line image has 
a point like emission (at NICMOS resolution) which is conci- 
dent with that of th e continuu m. The same happens for [Fe II] at 
1.6 urn ( Ma rconi et al . 2000). Regardless of these indications in 
order to avoid biases in the final results, we consider two cases, 
when the slits are centered on the emission line peak (coinci- 
dent with the continuum peak) and when they are not. 

The seeing of the observations can be estimated from 
the unresolved continuum source observed along the slit in 
the ISAAC observations. However, the unresolved continuum 
source in the J band is not as strong as in the K band and the 
estimate of the seeing can be uncertain by as much as 071. 
We parameterize the seeing as a gaussian with given FWHM 
(Full Width at Half Maximum) and we consider 3 cases for the 
FWHM values: 

- seeing 074 for PA1, PA2 and seeing 075 for PA3; 

- seeing 0745 for PA1, PA2 and seeing 0755 forPA3; 

- seeing 075 for PA1, PA2 and PA3; 

Inspection of Fig. [6] suggests that the ISBDs of Pa/8 and 
[Fe II] can be described with 4 components: a nuclear compo- 
nent, a disk component (the one associate d with the e xtended 
nuclear feature observed in Pao- bv lSchreier et al .119981) . and an 
extended one. The flattening and orientation of the disk com- 
ponents are fixed to q = 0.5 and 9 = 32.5, i.e. the values char- 
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Table 4. Fit results from the analysis of STIS data with i = 25 deg. 
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n Adopted ISBD. 

* Rescaled^ 2 with errors computed as AvJ 2 = Av, 2 + Avo 1 ■ 
c Systematic error adopted to renormalize^- 2 . 

d Average and rms of best fit parameter values. 

* Parameter was held fixed. 



acterizing the Paff feature, and the extended component is as- 
sumed circularly symmetric. We then parameterize the nuclear 
component according to the following 4 cases: 

1 . nuclear component circularly symmetric; 

2. nuclear component with PA = 32.5 deg and q = 0.5; 

3. nuclear component with PA = 122.5 deg and q = 0.5; 

4. 2 nuclear components with PA = 77.5, -12.5 deg and q = 
0.5. 

The adopted values of PA and q were chosen as described in 
the previous section. For each of the above cases we have used 
exponential (E), gaussian (G) and constant (C) functions. A 
surface brightness distribution is labeled by assigning a letter 
E, G or C according to the functional form used, a number 1-4 
(seeing case 1), 5-8 (seeing case 2) or 9-c (seeing case 3). The 
label is further characterized by 'b' if we are in the case when 
the emission line peak is not coincident with the slit centers. 
Therefore E5 indicates, for instance, exponential functions with 
surface brightness case 1 and seeing case 5 while G2b indicates 
gaussian functions in surface brightness case 2 and seeing case 
1. 



Sample results of the fitting procedure are shown in Fig. 1161 
for ISBD G7 and E7 and the best fit parameters for these ISBDs 
shown in Tab. [3] 

In Fig.^]we plot the isophotal contours of the ISBDs used 
in Figs. [2] and ^] These ISBDs have been chosen as rep- 
resentative because they provide the best fits of STIS veloci- 
ties (Tab. 0} and of ISAAC velocities and velocity dispersions 
(Tab. |6j. From the figure the centrally peaked components are 
clearly distinguished from the more extended ones. The central 
strong components of the ISBDs adopted for the STIS data is 
elongated along the jet axis, while the ISAAC ones are perpen- 
dicular to it. The more extended component in t he ISAAC data 
repres ents the disk like feature seen in Paff by ISchreier et all 
(1996) while it seems that this component is missing in [S III] 
which shows a different morphology from those of Pa/3 and 
[Fell]. 

4.3. Kinematical Fitting of the Rotation Curves 

In principle, in order to constrain the BH mass, one should 
compare the emission line profile predicted by the model 
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VIS1 VIS2 




-0.8 -0.4 0.0 0.4 -0.8 -0.4 0.0 0.4 

Pos. along slit [arcsec] Pos. along slit [arcsec] 



Fig. 18. Top panels: Fit of the STIS velocities along the slit. The empty squares with error bars represent the observed values 
while the solid lines connect corresponding model values. Dotted lines represent the contribution of the mass in stars to the 
rotation curve, i.e. what would be observed without a BH. 'G2' and 'E4' are the two ISBDs which provide the best fit models 
with the lowest^ 2 . Bottom panels: Velocity dispersions expected from the two best fit models compared with the observed ones. 
Notation as in the upper panels. 



with the obse r ved on e. However, as discussed at length by ments of the line profiles: the average velocity (v) and velocity 
iMarconi et al.l J2003ah . our approach is to compare the mo- dispersion <x defined as a 2 = (v 2 ) - (v) 2 . 
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Fig. 19. (a) x 1 contours for the joint variation of Mbh and M/L using the STIS data (ISBD G2). Contour levels are for^ 2 = 
X 2 ■ + 2.3, 4.61, 6. 17 corresponding to 68.3%, 90% and 95.4% confidence levels. The dot indicates the Mbh and M/L values for 
which the;^ 2 minimization was computed, (b) same as (a) but for Mbh - S. (c) Dependence ofx 2 on the adopted disk inclination. 
Dotted lines indicate the 68.3%, 90% and 95.4% confidence levels. 



Table 5. Effect of i variation on best fit parameters and mini- 
mum^ 2 for STIS data. All models adopt ISBD G2. Confidence 
levels 68.3%, 90% and 95.4% are at y 1 = x 2 mia + 10, 2.71, 4.0, 
for one interesting parameter ( Avni 1976). 



i log M BH log(M BH sin 2 9 xl: 
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7.68 


7.29 
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7.62 


7.32 
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7.56 
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163.7 


26.5 


55. 


7.52 
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165.0 
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7.55 


7.46 


164.7 


32.9 


75. 


7.64 


7.61 


161.9 


39.7 


85. 


8.17 


8.17 


160.9 


48.0 



This approach requires that the parameterized line profile, 
adopted in the fitting process, is able to reproduce t he shape of 
the ob served line profiles within the noise (see also Bart h et al.l 
200lJ. This condition is satisfied for STIS data were the ob- 
served line profile is well represented by a single gaussian, af- 
ter deblending the "blue" component which is obviously not 
emitted by circularly rotating gas. The above condition is also 
satisfied for the ISAAC data where the observed line profiles 
in the nuclear region are well matched with the parameteriza- 
tion described above. The Hermite expansion has the obvious 
advantage of allowing an excellent description of the line pro- 
files without introducing the instability and the larger number 
of parameters (6 instead of 5) typical of fits with two gaus- 



sian function. To our knowledge this is the first time that the 
Hermite expansion is applied to the analysis of emission lines. 

In order to be compared with the observations, model (v) 
and <x are computed taking into account the spatial resolution 
of the observations and the size of the apertures which are 
used to measure the observed quantities. The formulae used 
to compute velocities, widths and line surface brightness distri- 
butions and the details of their derivation, numerical computa- 
tion and m odel fitting are described in details in Appendix B of 
iMarconi e t al. (2003_a|)- 

In AppendixlAlwe describe a problem in the computation of 
kinematical quantities which might derive from subsampling. 
We also present a way to overcome this problem and therefore 
our computations are not affected by it. 

Thus, model (v) and cr depend on E, the ISBD of the emis- 
sion lines defined in the previous section, and on the following 
parameters: 

- xo,yo, the position of the kinematical center in the plane 
of the sky with the reference frame defined in the previous 
section; 

- i, the inclination of the rotating disk; 

- 0, the position angle of the disk line of nodes; 

- v sys , the systemic velocity of the disk; 

- M/L, the mass-to-light ratio of the stellar population; 

- M B h, the BH mass. 

Not all of these parameters can be independently determined 
by fitting the obser vations. It can be infer red from the equations 
in Appendix B of IMarconi et alJ J2003al) that M B h, M/L and i 
are directly coupled an d thi s has been discussed i n detail by 
iMacc hetto et all J 19971) and IMarconi et alJ J2003al) . To avoid 
problems deriving from this coupling we work here at fixed 
disk inclination i. 

4.3.1. STIS data 

We fit the STIS rotation curves using all the ISBDs deter- 
mined in Sec. 14. 2. H and, at first, we work with fixed inclination 
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Table 6. Fit results from the analysis of ISAAC data with i = 25 deg. 
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" Adopted ISBD. 

* Rescaled^ 2 with errors computed as AvJ 2 = Av, 2 + Avo 2 . 
r Systematic error adopted to renormalize^ 2 . 

d Average and rms of best fit parameter values. 
' Average and rms only on non-'b' ISBDs. 
f Average and rms only on 'b' ISBDs. 

* Parameter was held fixed. 
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Fig. 20. Top panels: Fit of the ISAAC velocities along the slit. The empty squares with error bars represent the observed values 
while the solid lines connects corresponding model values. Dotted lines represent the contribution of the mass in stars to the 
rotation curve, i.e. what would be observed without a BH. 'Gl' and 'E3' are the two ISBDs which provide the best fit models 
with the lowest^- 2 . Bottom panels: Velocity dispersions expected from the two best fit models compared with the observed ones. 
Same notation as in the upper panel. 



i = 25 deg. The STIS PSF used in the computations is obtained 
at A 9550A with TinyTIM jKrist & Hooklll999l) . The free pa- 
rameters of the fit (xQ,yo, logMBH, 0, v sys ) are obtained with 
a x 1 minimization. Since, as shown below, the contribution of 
the stellar mass to the gravitational potential is negligible we 
worked with fixed M/L=l in order to reduce the number of 
free parameters. We will show that this assumption is consis- 
tent with the data, when searching for the confidence limits on 
the parameters. During the fitting process, velocity corrections 
are applied to the VIS 1 and VIS2 data in order to take into ac- 
count errors on absolute wavelegth calibration. The best fit pa- 
rameters for all models with acceptable^ 2 are shown in Tab.0] 



and in Fig. ED (top panels) where we compare the observa- 
tions and the 2 models with the lowest^- 2 . In order to select the 
"good" models shown in Tab.|4]we first renormalize the y 2 fol- 
lowing a procedure al ready adopted bv lBarth et 

ID Good 

and 

Mar coni et alJ(l2003al) . The fit with the lowest^- 2 is provided by 
the model with ISBD "G2" which has x 2 kA = 1-27 (27.9/22). 
This ^ 2 ed value is already acceptable at the 95% confidence 
level given the number of degrees of freedom (22 dof) but we 
conservatively assume that there is an extra error in velocity, 
Avn, which has not been taken into account and that must be 
added in quadrature to the data. Avo might also represent an 
intrinsic spread in velocity due, for instance, to non-circular 
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local motions superimposed to the large scale rotation. Avo is 
chosen such the best fit model has X 2 e d~^- F° r model "G2", 
Avo = 10.4 km s~' and this is a small value compared to the 
errors on STIS velocities indicating the goodness of the fit. We 
then renormalize x 2 for a ll models using the above value of 
Avo and we consider acceptable only the models with a < 95% 
significance level, i.e. those with^ ed < 1.55 for 22 dof. For 
these models, the velocity shifts applied to the data are on aver- 
age -2.5kms _I (+2.6 rms) and 2. 6 km s" 1 (±2.7 rms) for VIS 1 
and VIS2 respectively. These are lower than the expected errors 
on th e absolute wavelength calibration (e.g. Kim Oui iano et alJ 
2003|. The models shown in Fig. ED are in good agreement 
with the observations as confirmed by the low found prior 
to rescaling. In particular the models nicely reproduce the steep 
velocity gradients and turnoffs observed in the rotation curves 
which are the signature of the BH. The velocity fall-off in VIS 1 
is slightly faster than keplerian but this disagreement appears in 
points with low signal-to-noise ratio. The best fit model is char- 
acterized by a BH mass of log(MBH/M ) = 8 and by a disk line 
of nodes with PA=161 deg. The dotted line in the figure shows 
the stellar contribution to the rotation curve, i.e. what would be 
observed without a BH. Clearly, the high velocity gradients can 
be explained only with the presence of a BH. The location of 
the kinematical center is -0 , .'05,-0705 with respect to the loca- 
tion of the slit centers. 

In order to estimate statistical errors associated with the 
best fit parameters, we present in Fig. I19h.b x 2 contours for 
the joint variation of Mbh - M/L and of Mbh - d. The^ 2 grids 
have been computed with ISBD "G2", fixing 2 interesting pa- 
rameters (Mbh, M/L or Mbh, ff) and varying the others to min- 
imize^ 2 . Before computing confidence levels we have renor- 
malized;^ 2 as described above so that the model with the lowest 
X 2 value has i y 2 ed =l. Confidenc e leve l s for 2 interesting param- 
eters are then found following lAvnil Jl976ll : x 2 - xt& + ^-3, 
4.61, 6.17 for confidence levels 68.3%, 90% and 95.4%, re- 
spectively. Statistical errors are +0.1 on log Mbh and ±6 on 6. 
Fig. 119b shows that it is not possible to determine M/L but only 
place an upper limit on it (M/L < 1.6). Also, the same figure 
justifies our choice of fixing M/L=l, which is within the lcr 
confidence contour. 

By comparing the results of the model fitting performed 
with different ISBDs we can have an estimate of the systematic 
errors associated with the choice of the ISBD. In the last row 
of Tab. |4]we show the average values of the fit parameters and 
the associated rms. It is highly significant that the systematic 
errors on all parameters and in particular on BH mass are very 
small when considering all the adopted ISB. The average BH 
mass is log Mbh = 7.88 with rms of 0.05. 

In conclusion, for i - 25 deg, the BH mass estimate from 
STIS data with ISBD "G2" is log Mbh = 8.0 ± 0.1 (±0.05 sys- 
tematic) and = 161 ± 6 (±6 systematic). 

In Fig. and Tab. we adopt ISBD G2 and we ex- 
plore how the minimum x 2 depends on disk inclination. As 
above, we have renormalized so that the model with the min- 
imum^ 2 value has^- 2 ed =l. The STIS data favour inclinations 
i < 48 deg at a 95% c onfidence level, a result already found by 
iMarconiet.al l fcOOll) . 



In Fig. 118b we plot the expected velocity dispersion for the 
2 best fit models plotted in Fig. II 8h . A constant intrinsic ve- 
locity dispersion of cr () = 160kms _1 was added to match the 
'plateau' observed at VIS1 (at 9500A the degradation of the 
spectral resolution introduced by the grating and slit width is 
~ 160kms _I , see Sec. l2.1> . In both cases an increase in cr is 
observed at the nucleus and this is roughly matched with un- 
resolved rotation only at VIS2. At VIS1 unresolved rotation 
explains only ~ 80% of observed cr nuclear increase but the 
adopted cr = 160kms _1 is smaller than the contribution from 
unresolved rotation. 

4.3.2. ISAAC data 

We fit the ISAAC rotation curves using all ISBDs determined 
in Sec. 14.2.21 We proceed as in the previous section and we 
work at fixed inclination, i = 25 deg. The PSF used in the com- 
putations is a gaussian with FWHM equal to the adopted see- 
ing. The free parameters of the fit (xo,yo, log Mbh, log M/L, 
9, v sys ) are obtained with a x 2 minimization. During the fit- 
ting process, velocity corrections are applied to the PA1, PA2 
and PA3 data in order to take into account errors on absolute 
wavelegth calibration. The best fit parameters for all models 
with acceptable^ 2 are shown in the first part of Tab.[6]("Fit of 
velocity") and in Fig.|20|we compare the observations and the 
2 models with the lowest^ 2 . 

In order to select the "good" models shown in Tab. [6] we 
renormalize the x 2 following the approach described in the 
previous section. The fit with the lowest x 2 is provided by 
the model with ISBD "Gl" which has ^ 2 ed =2.62 (91.9/35) and 
must be then renormalized with Avo = 6.3 km s _1 . We then con- 
sider as acceptable all the models with a significance of 95%, 
i.e. those with^ 2 ed < 1.37 (35 dof). For the models in Tab. [6] 
the average velocity shifts applied to the data are -1.0±3.4, 
-18.7±4.4 and 10.9±5.7 for PA1, PA2 and PA3 respectively. 
These are in very good agreement with the expected errors on 
the absolute wavelength calibration (Marc oni et al.l200ll) . 

As for the STIS data, the models reproduce the observed 
velocities and, in particular, they match the steep velocity gra- 
dients and turnoffs observed in the rotation curves which are 
the signature of any BH. The best fit model is characterized 
by a BH mass of log(M B H/M ) = 8.14 and by a disk line of 
nodes with PA=164 deg. The location of the kinematical center 
is -0.08,-0. 16 with respect to the location of the slit centers. 

In the lower panel of Fig.|20|we plot the expected veloc- 
ity dispersion for the 2 lowest x 2 models. A constant intrinsic 
velocity dispersion of cr = 70kms _I was added to match the 
'plateau' observed at PA1. In both cases an increase in cr is 
observed at the nucleus but this is clearly not matched by the 
models with unresolved rotation. 

In order to verify if the observed cr increase can be repro- 
duced by the models we also perfomed the fits by consider- 
ing velocity and velocity dispersion at the same time. Best fits 
acceptable to the 95% level (i.e., with 75 d.o.f., renormalized 
X 2 ed < 1 .28) are shown in the second part of Tab.|6]("Fit of ve- 
locity and velocity dispersion") and in Fig.|^we compare the 
observations and the 2 models with the lowest^ 2 . Though the 
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Fig. 21. Fit of the ISAAC velocity and velocity dispersion along the slit. The deviant velocity dispersions in PA1 and PA2 (points 
not marked with an empty square) were not included in the fit. Notation is the same as in the previous figure. 



fit of velocity is not as good as before, the observed curves are 
still acceptably reproduced by models. However this time mod- 
els are also able to match the observed sigma increase showing 
that it can be explained by unresolved rotation. The reason for 
this result is that in the previous case the center of rotation was 
pushed far from the surface brightness peak in order to better 
reproduce velocities. Now the center of rotation is almost co- 
incident (within ~ (X'01) with the line emission peak (note that 
in the case of ISBDs 'b' the emission line peak is offset from 
the slit centers by approximately 0.1, -0.02 and this is where 
the kinematical center is). 

In Fig. [22] we show the y 2 contours for joint variation of 
Mbu,M/L and Mbh,& in order to estimate the statistical errors. 
The models have been computed adopting ISBD G7, fixing 2 



interesting parameters and varying the others. The^ 2 has been 
renormalized as before. 

The statistical errors associated to the BH mass determina- 
tion are +0.05 in log Mbh comparable to the systematic errors 
due to the adopted ISBD. As for the fit of the STIS data Fig.l22l 
shows that it is not possible to determine M/L but one can only 
place an upper limit of M/L < 1 . The statistical errors associ- 
ated to 9 are +3 deg. 

By comparing the results of the model fitting performed 
with ISBDs we can have an estimate of the systematic errors 
associated with the choice of the ISBD. In Tab. [6] we also 
show the average values of the fit parameters and the associ- 
ated rms. It is highly significant that the systematic errors on 
all parameters and in particular on BH mass are very small. 
But the most important thing is that whether one fits veloc- 



24 



A. Marconi et al.: The Black Hole in Centaurus A 




log M B] , [MJ log M BH [M .] Inclination [deg] 



Fig. 22. (a) x 1 contours for the joint variation of Mbh and M/L using the ISAAC data and fitting velocities and central velocity 
dispersions (ISBD G7). Contour levels are for^f 2 - x 2 ■ + 2.3, 4.61, 6.17 corresponding to 68.3%, 90% and 95.4% confidence 
levels. The dot indicates the Mbh and M/L values for which the^ 2 minimization was computed, (b) same as (a) but for Mbh - 8. 
(c) Dependence of^- 2 on the adopted disk inclination. Dotted lines indicate the 68.3%, 90% and 95.4% confidence levels. 



Table 7. Effect of i variation on best fit parameters and 
minimum x 2 f° r ISAAC data. All models adopt ISBD G7. 
Confidence levels 68.3%, 90% and 95.4 % are at y 1 = x 2 oiD + 1 0, 
2.71, 4.0, for one interesting parameter ( Avni 1 9761) . y 2 's have 
been rescaled as described in the text with AVo = 4.6 km s" 1 . 
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ity or velocity and sigma, the BH mass is the same. Indeed 
log M B h = 8.14 ± 0.04 (average over different ISBDs, velocity 
only), log M B h = 8.11 + 0.08 (average over different ISBDs, 
velocity and velocity dispersion). Also all fl's are consistent. 

In conclusion, for i = 25 deg and ISBD G7, the BH mass 
estimate is logM B H = 8.03 + 0.05 (+0.06 systematic) and 8 - 
165 ± 3 (+1 systematic). 

In Fig. 122b and Tab. [7] we explore how the minimum x 2 
varies with varying disk inclination for the models adopting 
ISBD G7. As for the STIS case, we have renormalized x 2 ' & 
so that the model with the minimum^ 2 value has^ 2 ed =l. The 
ISAAC data clearly favour inclinations i < 35 deg at the 95% 
confidence level. This limit on i is smaller than the one derived 



from STIS data however, given the higher SNR of the ISAAC 
data and the larger number of data points, we take it as the 
upper limit on i. 

5. Discussion 

In the previous sections we have described in detail the avail- 
able datasets and the application of the gas kinematical method 
which provides consistent results from independent analysis of 
STIS and ISAAC data. In this section we analyze the influence 
of the ISBD on the final results and we discuss in detail the 
results obtained using STIS and ISAAC data. We then verify 
the crucial assumption that the gas is in a thin, circularly rotat- 
ing disk by comparing observed and model velocity dispersion. 
We caution about incorrect computations of the model kine- 
matical quantities which might be affected by a subsampling 
problem (see also AppendixlAli or by adopted ISBD which are 
too smooth. We also compare observed line profiles, /13 and 
parameters with model values. Finally, after comparing the gas 
kinematical estim ate of Mbh with the stellar dynamical one by 
Silg e et all J2005), we discuss Mbh in the framework of the 
correlations between BH mass and host galaxy structural pa- 
rameters and we conclude this section by estimating the maxi- 
mum size of a cluster of dark stellar remnants, the most likely 
alternative to a supermassive BH. 

5.1. Influence of the intrinsic surface brightness 
distribution 

One of the worries in gas kinematical measurements of BH 
masses has always been the role of the intrinsic emission 
line surface brightness distribution, since it is the weight 
of the velocity field in the averaging over apertures (e.g . 
iMacchetto et al.ll997tlBarth et al.l200UlMarconi et al.l2003 d) . 
In the previous section we have clearly shown that the adopted 
ISBD, provided that it reproduces the observed one within the 
errors, does not affect the final BH mass estimate. Indeed sys- 
tematic errors on M B h due to the adopted ISBD are of the order 
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of 0.05 and 0.08 in log M BH for STIS and ISAAC data respec- 
tively. However it should be remembered that the ISBD has an 
important effect on the quality of the velocity fit. 

We can now compare the res ults obtained here with those 
obtained bv lMarconi et al.l 112001 ) who adopted a simple con- 
stant ISBD and a double exponential. From the ISAAC data 
we obtain logM B n = 8.11 + 0.08 (average and rms over all 
ISBDs ) with disk incl ination i = 25. With the same disk incli- 
nation, Marc oni etaD d200ll) obtained logM BH = 8.31 + 0.1 
(constant ISBD) and logM BH = 8.25 ± 0.1 (double expo- 
nential ISBD). The data are consistent within the errors. The 
slightly small er values f ound here are explained with the fact 
that Mar coni et al.l J200ll) did not consider the stellar mass 
in their fitting. Rega rding the PA of the disk line of nodes 
iMarconi etaiTll200ll) found -15+5 while we find 164 ± 2, i.e. 
-16 + 2. Therefore our new analysis impr oves the accurac y of 
the BH mass estimate with respect to lMarconi et alJ |2001 ) but 
does not substantially change their results. 
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5.2. Comparison of results from fitting STIS and 
ISAAC data 

In the previous sections we presented the results of fitting sepa- 
rately STIS and ISAAC data. Since the fits are completely inde- 
pendent and data have spatial resolutions differing by an order 
of magnitude, the comparison of the results provides an im- 
portant check of the robustness of the gas kinematical method 
in measuring BH masses. Moreover, the comparison will show 
whether it is possible to obtain an accurate estimate of the BH 
mass even with lower S/n data like those obtained with STIS. 

We did not perform a joint fit of the ISAAC and STIS data 
because this would be strongly biased toward the ISAAC data 
since the number of points is larger and their relative errors are 
smaller. 

The best fit of STIS data with ISBD G2 (see Tab. g) 
provides log(M BH /M ) = 8.00 ± 0.11 (statistical) +0.05 
(systematic). This is perfectly in agreement with the best 
fit result of ISAAC data where, from ISBD G7 (Tab. |6}, 
log(M BH /M ) = 8.03 + 0.05 (statistical) +0.08 (systematic). 
The agreement is still good when comparing the best fit val- 
ues of the BH mass averaged over all the adopted ISBDs. STIS 
data provide log(M BH /M G ) = 7.88 ± 0.05 while from ISAAC 
log(M BH /M Q ) = 8.11 + 0.08. The log ratio of the BH masses 
is 0.2 + 0.1 which is consistent with at the 2cr level. 

Apart for the BH mass, the separate fits recover consistent 
estimates of the position angle of the disk line of nodes. From 
STIS (G2), the best fit provides 6 = 161 +6 (statistical) ±6 
(systematic) while from ISAAC (G7) it results 6 = 168 + 2 
(statistical) ±2 (systematic). The agreement is even improved 
when comparing the best fit values of 9 averaged over all the 
adopted ISBDs: = 168 ± 6 (STIS) and 6 = 164 + 2 (ISAAC). 

In Fig. [23] we compare the positions of the slit centers and 
emission line peaks with respect to the kinematical center posi- 
tion derived in the fitting. Since in each fit we leave the coordi- 
nates of the kinematical center free to vary, in order to perform 
the comparison we consider a reference system in which the 
kinematical center is at (0,0). In figure, slit centers are indi- 



Fig. 23. Positions of slit centers (empty symbols) and emission 
line peaks (filled symbols) compared to the kinematical center 
located at Aa = 0, A6 = 0. Ellipses indicate the lcr errors. The 
dashed line is the PA of the jet observed at radio and X-ray 
wavelengths. 

cated by empty big symbols while the corresponding positions 
of the emission line peaks are indicated by smaller filled sym- 
bols with the same shape. In the STIS case, the slits were po- 
sitioned blindly on the nucleus position which was computed 
from NICMOS K images comparing the location on the con- 
tinuum peak with respect to the acquisition star. The ISAAC 
slits were directly centered on the location of the K continuum 
peak. During the fitting of the ISAAC data we considered two 
cases, one where the emission line peak was coincident with 
the slit centers (i.e. with the K continuum peak as suggested by 
the NICMOS Paa image), and the other, labeled 'b', where the 
line emission peak was not coincident with the slit centers. In 
all cases slit centers are closer than 0. 1 arcsec from the kine- 
matical center. But the most interesting indication comes from 
the location of the emission line peaks which are, within the 
errors, located along the jet axis, a result not obvious a-priori. 
This is consistent with two possible pictures: in the first one 
the jet excites line emission while in the second one the jet rep- 
resents the symmetry axis of th e ionization cone, according to 
the un ified model of AGNs (e.g. lAxon et alJl998tlCapetfi et all 
1999). 

5.3. Is the gas circularly rotating? Velocity dispersion 
and line profiles tests. 

The power of gas kinematical measurements is their conceptual 
simplicity compared to the complex analysis required by stellar 
dynamical modeling for which the re are serious issue regardin g 
indeterminacy dValluri et alJ2004lCretton & Emsellerrj20 04). 
However, gas kinematical measurements are based on the as- 
sumption of Keplerian disks, i.e. thin disks circularly rotating 
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under the influence of the gravitational potential due to stars 
and BH. In some cases the gas is clearly not rotating circularly, 
for instance when the motions are affected by jet -cloud inter- 
actions (e.g. lAxon et al.lfl998t ICapetti et alJll999h but in other 
cases the question is more subtle. 

The comparison between observed and model velocity dis- 
persion has become one of the most used tests on the relia- 
bility of the circular rotation assumption: velocity dispersions 
larger than expected from unresolved rotation could be an in- 
dication of non circular motions which could invalidate the BH 
mass estimate. Many authors have found that the observed ve- 
locities could be well fit with Keplerian disk models but the 
gas velocity dispersion was larger, and in some cases much 
larger than expected from unresolved rotation. This has caused 
big concerns on the reliability of the gas kinematical m ethod 
(e.g. Ivan der Marel & van den Boschlll998l iBarth et alJhoOlt 
Verd oes Kleiin et al.l2000H2002HCappellari et al.l2002l) . 

Before focussing on Centaurus A, we show that these con- 
cerns might have been excessive since a small model velocity 
dispersion could easily result from (i) a non-accurate compu- 
tation or, more likely, from (ii) an adopted ISBD which is too 
smooth. 

In AppendixlAlwe show the effect of a coarse sampling in 
the computation of the model kinematical quantities. This ef- 
fect is particularly important when the ISBD varies on spatial 
scales which are much smaller than the PSF sizes and Fig. IA.ll 
shows how one could easily underestimate the model velocity 
dispersion by a factor ^ 1.5. This effect is also very subtle be- 
cause it is not manifested in the model computed line fluxes 
and velocities. We thus caution about the sampling chosen in 
computations and in Appendix|A]we present a possible way to 
overcome this potential problem. 

We now focus on the effects of different ISBs which provide 
the same light profile along the slit after being convolved with 
the instrumental response. We consider three different ISBs la- 
beled A, B and C. All of them are circularly symmetric on the 
plane of the sky and are described by a sum of two exponentials 
as 

/(r) = / e-'- / '-°+/ 1 e-'- / '-' (9) 

ISBD A has h/h = 0.01, r = 0.05", and n = 1.00". B 
has Ii/Io = 0.0067, ro = 0.08" but the second exponential 
is substituted with a constant. Finally C has h/Io = 0.0008, 
ro = 0.01", and r\ = 0.69". The 'observed' system is like 
Centaurus A but with M B h=1O 8 M , M/L-l, and the gas disk 
is inclined by 45 deg with respect to the line of sight. The rotat- 
ing gas disk, which has an intrinsic constant velocity dispersion 
of 70 km s , is observed with ISAAC with 0.5" seeing. The slit 
is placed along the major axis of the disk. In Fig. [23] we show 
the intrument convolved surface brigtness, velocity and veloc- 
ity dispersion along the slit. As noted in previous sections, it 
is clear that the observed surface brightness and velocity along 
the slit (and consequently Mbh) are little sensitive to the choice 
of the ISBD. The reason for this weak dependence of the instru- 
ment convolved velocity on the ISBD is that, close to the BH 
location where the weighting by the ISBD is larger, rotation 
is not resolved and the contribution from positive and nega- 
tive velocities is canceled. However, the effects of the ISBD 



are readily apparent on the gas velocity dispersion. Suppose 
that in the analysis of the data one uses ISBD A and finds that, 
while the observed velocities are well reproduced by the model, 
the gas velocity dispersion is underestimated by a factor =2 2. 
One might then be tempted to conclude that the gas has a high 
intrinsic velocity dispersion. However, using ISBD C (or even 
a less extreme one) one would reach the opposite conclusion 
regarding the validity of the Keplerian disk assumption. 

It is not possible to draw general conclusions here because 
each case must be analyzed separately bearing in mind the pos- 
sibility that the adopted ISBD does not reflect the real one. 
However we caution against the possibility that many of the 
previous gas kinematical analyses finding inconsistent gas ve- 
locity dispersions might be affected either by a subsmapling 
problem or by the adoption of a too extended ISBD. 

We now focus on Centaurus A where the STIS and ISAAC 
datasets are well suited to this discussion. In the case of ISAAC 
data, the model fitting of velocities does not reproduce the ve- 
locity dispersion and the reason is that, in order to find the best 
match of the rotation curves, the kinematical center is pushed 
away from the emission line peak. However, the joint fit of ve- 
locity and velocity dispersion shows that the circularly rotating 
model can well account for both velocity and velocity disper- 
sion at the same time (note that the Mbh estimate is not changed 
by this). 

The observed velocity dispersion in the STIS datasets is 
larger than expected from pure rotation but the previous discus- 
sion on the effects of the ISBD on the model velocity dispersion 
might suggest that it is possible to match the observed STIS 
velocity dispersion with an appropriate choice of the ISBD. 
However the main ISBD component used for STIS is not lo- 
cated on the BH but 0V 1 away, therefore changing its scale 
radius would not produce any appreciable effect on final model 
values. Instead, one should place a bright and very peaked 
ISBD component on the BH location. In Fig. [25] we compare 
model kinematical quantities obtained using G2 and two addi- 
tional ISBDs obtained by adding an exponential component to 
G2 at the location of the BH. These exponential components 
have h = 5000, r = 0.02" and I = 20000, r = 0.01" re- 
spectively (see also Tab. [2] for the parameters describing G2). 
Surface brightnesses and velocities are very little affected by 
this additional ISBD component located at the position of the 
BH. However the gas velocity dispersion is increased and it is 
now possible to reproduce the peak value of the observed veloc- 
ity dispersion. The observed cr profile along the slit is slightly 
wider than the model one nonetheless, given the signal-to-noise 
ratio of the data, the agreement between models and observa- 
tions can now be considered satisfactory. The intrinsic constant 
velocity dispersion required by the data, cr = 160 km s -1 , is 
smaller than that expected from circular rotation and rotation 
around a BH is needed to explain the increase of <x in the nu- 
clear region. The measured Mbh is perfectly consistent with 
the results from the ISAAC dataset. This is a clear indication 
that even the presence of a constant intrinsic velocity dispersion 
(cr = 160kms _1 ) does not necessarily invalidate the Mbh es- 
timate. In conclusion it is important to emphasize that a large 
velocity dispersion rise in the nucleus of a galaxy cannot be a 
manifestation of the ISBD in a pure stellar model. Thus their 
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Fig. 24. Instrument convolved kinematical quantities computed with different ISBDs A,B and C (see text for explanation). The 
ISBDs have been chosen to provide very similar observed surface brightnesses along the slit. 




3 



E I i i i i I i i i i I i , i i I i a 



400 



300 



200 



100 



400 



300 



200 



100 



-0.8 -0.4 0.0 



-0.8 -0.4 0.0 
Position along the slit [arcsec] 




-0.8 -0.4 0.0 



Fig. 25. Comparison between observed kinematical quantities from STIS (points with error bars) and model computations. From 
left to right panels represent surface brightness, velocity and velocity dispersion along the slit. The top and bottom panels are 
for VIS1 and VIS2, respectively. Models were computed using G2 and two additional ISBD contructed by adding to G2 an 
exponential component with ro = 0.02" and ro = 0.01" at the position of the BH (see text for more details). 



very existence is prima facie evidence for the presence of a BH 
provided the gas motions are dominated by rotation. However, 
due to the sensitivity of the observed sigma to the ISBD, mod- 
eling a sigma peak is model dependent and therefore results in 
weak constraints on the precise mass of the BH. 

In Fig. |26]we compare the model line profiles computed 
with ISBD G7 with the observed ones with ISAAC. Such a 
comparison is meaningless with STIS because of the much 
lower spectral resolution. The match between observed and 



model line profiles is good given that the fitting was performed 
on velocity and velocity dispersion and not on the line pro- 
file itself. G7 is the ISBD which provides the best fit of ve- 
locity and velocity dispersion. We have also searched for the 
ISBD which provides the best match of the line profiles among 
those in Tab.|6] A% 2 analysis indicates that this is given by E5 
and indeed also a visual inspection of the figure confirms this 
result. Model line profiles, which are only affected by unre- 
solved rotation, are able to reproduce, at least qualitatively, the 
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Fig. 26. Top panels: observed ISAAC line profiles are compared with model ones computed with the best ISBD G7. Numbers 
on the top right corners indicate the row where the spectrum was extracted (175 is the location of the slit center). Model profiles 
have been rescaled in surface brightness to have the best match with observed ones. Bottom: as in top panels, but for ISBD E5 
which provides the best matching model line profiles. 
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Fig. 27. Top panels: observed Hermite parameters and h\ compared with model ones computed with the best ISBD G7. 
Numbers on the top right corners indicate the row where the spectrum was extracted (175 is the location of the slit center). 
Bottom: as in top panels, but for ISBD E8 which provides the best match to the observations. 
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Table 8. Comparison of BH mass measurements with gas kine- 
matics and stellar dynamics 



i 


Gas Kinematics" 
M BH ;t 2 /dofO 


Stellar Dynamics 
M B h X 2 /dof 


20 
45 
90* 


1.9^ x 10 8 
4.8 + ^ x 10 7 
x 10 8 


41/41 
44/41 
130/41 


1.5 + ^xl0 8 ... 
1.8^ xlO 8 ... 
2.4** x 10 s ... 



" Rescaled,^ 2 from tabled 



For gas kinematics i = 85 deg. 

features which are observed in the data li ke asymmetrie s and 
humps (see also the discussion in iMacchetto et al.1l 19971) . It is 
thus clear that line profiles are consistent with the rotating disk 
model and that it is possible to discriminate among different 
ISBD on the basis of the line profiles. 

In this paper we have, for the first time, introduced the 
Hermite expansion in the analysis of emission line profiles. It is 
therefore natural to veirfy whether the model values of /13 and 
/14 are consistent with observed ones. In Fig. \^\ we perform 
this comparison using the best fitting ISBD G7. The agreement 
is not very good but it must be kept in mind that the fit was 
not optimized to reproduce hj and /14. Qualitatively the model 
reproduces the shapes of the h.3 and h\ variations along the slit, 
though it does not reproduce the amplitudes. As in the previous 
case, we have searched for the best fit model which provides the 
best match to and I14 and this results to be that generated with 
ISBD E8. The match is much better than previously though not 
yet as good as the match of the observed velocity and velocity 
dispersion. As in the previous case we can conclude that the 
comparison with hj, and hi, can be used to discriminate among 
different ISBDs but this is beyond the scope of this paper. 

5.4. Comparison with Stellar Dynamical 
Measurements 

ISilge et alJ fe005t very recently presented a stellar dynamical 
measurement of the BH mass in Centaurus A. They use ax- 
isymmetric 3-integral stellar dynamical models based on the 
orbit superposition method by Schwarzschild. They estimate a 
BH mass of 2.4^ 2 x 1O 8 M for edge-on models (i.e. where 
the principal plane of the potential is edge-on, i = 90 deg, with 
respect to the line of sight), l.S^ x 10 8 M o for i = 45 deg, and 
1.5^2 x 10 8 M Q for i = 20 deg. In the Keplerian disk approx- 
imation, the gas disk is rotating in the galaxy principal plane. 
Therefore the inclination of the principal plane is directly the 
inclination of the gas disk. Thus, in table [5] we compare gas 
kinematics and stellar dynamical BH mass measurements com- 
puted with the same inclination. 

It is readily apparent that there is an excellent agreement 
between the two measurements if we consider the i = 20 deg 
case while, at intermediate inclination i = 45 deg, the measure- 
ments are discrepant by at least a factor of 2. The good agree- 



ment for i = 90 deg is not significant because the^ 2 in the gas 
kinematical analysis is a factor of 3 worse than at other incli- 
nations. As discussed in the previous sections, gas kinematics 
favour low inclinations of the gas disk and, indeed, the model 
at i = 20 deg has the lowest x 2 among those in the table. On 
the contrary, stellar dynamical models suggest an edge on prin- 
cipal plane which is completely excluded by gas kinematical 
models. 

The stellar dynamical measurements are potentially af- 
fected by severa l problems, first of all, as discussed in 
ISilee et all J2005I) . the possibility that Centaurus A is not 
axisymmetric. Centaurus A has a spherical simmetry in the 
nuclear region but it seems moderately triaxial in the outer 
regions. At variance with gas kinematics which can work 
only with the rotation curves from the nuclear region, stel- 
lar dynamics needs extended data to constrain st ellar orbits 
(e.g. lValluri et alJl2004l: ICretton & Emsellerrl2004l) . As a con- 
sequence, Sil ge et al.l d2005b conclude that it is unlikely that 
they can constrain the actual inclination of Centaurus A using 
an axisymmetric modeling procedure. 

Moreover, the inclination of the gas disk is not necessarily 
the same as that of the galaxy principal plane. In the region 
where we obtain the gas kinematics the gravitational potential 
is completely dominated by the BH and the gas might lie in a 
plane which is not coincident with the principal plane of the 
galaxy. 

The above considerations allow us to relax the constraint 
on having the same inclination for gas and stellar kinematical 
models and we can thus conclude that the two methods provide 
results in excellent agreement. This is the first time, for an ex- 
ternal galaxy, that there are consistent BH mass measurements 
with gas and stellar dynamics. 

This agreement should not be surprising since it is found 
even for the closest supermassive BH, that in the Galactic 
Center. Notwithstanding the comp lex environment, g as kine- 
matical measurements from Genzel & Townes ( 1987) provide 
a Mbh measurement which is i n excellent agreemen t with that 
of the stars (see, e.g., Fig. 11 of lSchodel et al. 2003). The only 
problem of the gas is that it is not present in the inner 0. 1 pc 
and therefore the constraints it poses on the size of the massive 
dark object are much weaker that those of the stars which can 
probe the gravitational potential more than 2 orders of magni- 
tude closer to the location of Sgr A*. 

5.5. Correlations with galaxy properties 

The BH mass estimate in Centaurus A can be summarized as 
Mbh = (1.1 ± 0.1) x 10 8 M o for an assumed disk inclination 
of i = 25 deg or M BH = (6.5 + 0.7) x 1O 7 M for i = 35 deg. 
We now compare these two values with the known correlations 
with host spheroid lum inosity and mass (we consider for the 
analysis the results by Marconi & Hunt 2003) and with the stel- 
lar ve locity dispersion (Tre maineeta?T2002t[Ferrarese & Ford 
2005). 

We use the K band Mbh - i-sph correlation by 
iMarconi & Hunt! J2003I) who consider only galaxies with "se- 
cure" Mbh determination. The K band total luminosity of 
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Fig. 28. Maximum core radius from STIS observations (ISBD 
G2 and i = 25). 

Centaurus A is estimated bv lMarconi & Hunt! J2003I) as M K = 
-24.5 (L K = 1.3 x 10 11 L 0>K ) and the BH mass expected from 
the correlation is ~ 2.8 x 1O 8 M , in agreement with our results 
if one takes into account observational errors and the intrinsic 
scatter of the M BH - £ S nh correlation which is 0.3 in log M sli (a 
factor of 2). |m arconi &Hunt (2003) also found that the virial 
mass [M sp h = 3R e crl/G] is correlated with Mbh. Using their es- 
timate of M sp h for Centaurus A (M sph = (5.6 ± 1.5) x 1O 1O M ) 
and the correlation which considers only galaxies with "secure" 
Mbh determination, the expected BH mass for Centaurus A is 
~ 1 .4 x 1O 8 M , again in excellent agreement with our measure- 
ment. Therefore we can conclude that Mbh in Centaurus A is 
in very good agreement with L sp h and M sp h. 

The luminosity weighted stellar velocity dispersion of 
Centaurus A using an a perture of 60" par allel to the dust lane 
is a = 138 + lOkms -1 (Sil ee et alJl2005l) . Assuming that this 
is a good approximation for cr e (the same quantity but aver- 
aged in an aperture R e , the galaxy effective radius), we can 
estimate the expected BH mass according to the M B h - cr e 
correlation. Adopting the correlation parameters estimated by 
iTremaine etafl J2002I) the expected BH mass is ~ 3.0 x 1O 7 M 
while with the version of the correlation by Fe rrarese & Ford 
J2005I) M BH ~ 3.0 x 1O 7 M . The BH in Centaurus A is thus, 
apparently, a factor 2-4 more massive than expected from the 
correlations. However, if the scatter of the intrinsic scatter of 
the Mbh - o"e correlation is taken into account (0.3 in log Mbh), 
then the discrepancy can be fully accounted for. In conclusion, 
if the Mbh _ o"e is not a perfect correlation, like it is reasonable, 
then the BH in Centaurus A is not a deviant point. 

5.6. The Size of the Dark Cluster 

The kinematical analysis which leads to the discovery of a BH 
and measurement of its mass in reality produces evidence only 



Fig. 29. Loci of dark cluster lifetimes (yr) in the mass-density 
plane (after lMaoJ 19981) . 



for the existence of a massive dark object. The size of this ob- 
ject is usually constrained to be smaller that the spatial resolu- 
tion of the observations but in most cases , except for our own 
Galactic Center (e.g. Scho del et alJl2003l) . this is not enough 
to exclude that the massive dark object is not a cluster of dark 
objects. 

In the present paper we can take advantage of the high 
spatial resolution achieved with HST to estimate a tight up- 
per limit to the size of the putative dark cluster. We assume 
that the dark mass is distributed according to a Plummer model 
(e.g. iBinnev & Tremainelll987h and we explore the x 2 varia- 
tion for fixed values of the core radius, R C0K , minimizing all 
the other free parameters of the fit xq, y , log M core (core mass), 
9, v sys . Results, using ISBD G2 and i = 25deg, are shown in 
Fig. [28] As before, \ 2 is rescaled in order to have = 1 

for the best fitting model and confidence levels are computed 

I 1 I I 

following Avni ( 1976). To the 90% confidence level, the core 
radius is less than 07036 (0.6 pc). This value increases to 0704 
(0.7 pc) for a disk inclination of i = 35 deg. This implies that 
the putative massive dark cluster must have a minimum den- 
sity ofpBH = M BH /(4/37r^ 01 . e ) = 1.1 x 10 8 M o pc -3 . This value 
is almost a factor 6 larger than that found for M87 which, af- 
ter our o wn galaxy an d NGC 4258, was among the best cases 
for a BH. Maoz ( 1998) computed the lifetimes of massive dark 
clusters with given mass and density before their collapse to a 
supermassive BH. This lifetime, if much shorter than the age 
of the universe, indicates that the most likel y alternative i s that 
of a BH. Fig.|29lis an update of Fig. 1 bylMaoTl il 19981) with 
the compilation by Marconi & Hunt ( 2003). The stars mark the 
new location of Centaurus A based on HST STIS observations 
for inclination i = 25 and 35 deg. For comparison, we also 
show the previous location of Centaurus A based on ground 
based ISAAC observations only. The solid lines are the loci of 
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constant cluster lifetimes computed bv lMaozl ( fl9 98). The life- 
time of the putative dark cluster in Centaurus A is still larger 
than 10 11 yr and comparable to the age of the universe indi- 
cating that a cluster of dark objects cannot be ruled out. Still, 
Centaurus A ranks now among the best cases for a BH indicat- 
ing that much work is still waiting to secure BH mass deter- 
minations and the proof on the existence of supermassive BHs. 
To improve in this direction one needs spatial resolution of the 
order of 1 mas such as those which can be rea ched with new 
genera tion interferometers like the VLTI (e.g. iMarconi et alJ 
l2003bl) . 

6. Summary and Conclusions 

We have presented new HST Space Telescope Imaging 
Spectrograph observations of the nearby radio galaxy NGC 
5128 (Centaurus A). Since the galaxy nucleus is reddened by 
Ay > 7 mag, we used the the redder emission line accessible 
from HST, [S III]/19533A, to study the kinematics of the ion- 
ized gas at 07 1 spatial resolution. 

The STIS data were analized in conjunction with the 
ground-bas ed near-infrared Very Large Telescope ISAAC spec- 
tra used by Marconi et alJ d200lh to infer the presence of a su- 
permassive black hole and measure its mass. The ISAAC data 
were re-analyzed and we introduced the use of the Hermite 
expansion for the kinematical analysis of emission lines. The 
Hermite expansion allows an excellent fit of the observed line 
profiles with a smaller number of parameters and more stability 
than two gaussian functions. 

The two sets of data have spatial resolutions differing by 
almost a factor of five but provide independent and consistent 
measures of the BH mass, which are in agreement with our pre- 
vious estimate based on the ISAAC data alone. The gas kine- 
matical analysis provides a mass ofMeH = (1.1 ±0.1)x 1O 8 M 
for an assumed disk inclination of i = 25 deg or Mbh = 
(6.5 + 0.7) x 1O 7 M for i = 35 deg. Our new analysis im- 
proves the accuracy of the BH mass estimate with respect to 
Mar coni et al.l d200ll) but does not substantially change their 
results. 

We performed a detailed analysis of the effects on Mbh of 
the intrinsic surface brightness distribution of the emission line, 
a crucial ingredient in the gas kinematical analysis. We estimate 
that the associated systematic errors are no larger than 0.08 in 
log Mbh, comparable with statistical errors and indicating that 
the method is robust. However, the adopted intrinsic surface 
brightness distribution has a large impact on the value of the 
gas velocity dispersion. A mismatch between the observed and 
model velocity is not necessarily an indication of non-circular 
motions or kinematically hot gas, but is as easily due to an in- 
accurate computation arising from too course a model grid, or 
the adoption of an intrinsic brightness distribution which is too 
smooth. 

The velocity dispersion in the ISAAC spectra is matched 
with a circularly rotating disk and also the observed line pro- 
files and the higher order moments in the Hermite expansion 
of the line profiles, /13 and /14, are consistent with emission 
from such a disk. The velocity dispersion in the STIS data is 
slightly larger than expected from rotation however this mis- 



match can be reconciled when taking into account the signal-to- 
noise ratio of the data and the possibility of using more peaked 
intrinsic surface brightness distributions described in the text. 
Nonetheless, STIS data provide an estimate of the BH mass 
consistent with ISAAC indicating that even an intrinsic veloc- 
ity dispersion of the gas does not invalidate gas kinematical BH 
estimates. 

To our knowledge, Centaurus A is the first external galaxy 
for which reliable BH mass measurements from gas and stel- 
lar dynamics are available. The agreement between our Mbh 
gas kinematical estimate and the similar estimate from stellar 
dynamics thus strengthens the reliability of both methods. 

The BH mass in Centaurus A is in excellent agreement with 
the correlation with infrared luminosity and mass of the host 
spheroid but is a factor ~ 2 - 4 above the one with the stellar 
velocity dispersion. But this disagreement is not large if one 
takes into account the intrinsic scatter of the Mbh ~ o"e correla- 
tion. 

Finally, the high HST spatial resolution allows us to con- 
strain the size of any cluster of dark objects alternative to a BH 
to r. < O/035 (- 0.6 pc). Thus Centaurus A ranks among the 
best cases for supermassive BHs in galactic nuclei. 
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Appendix A: Computation of model kinematical 
quantities: the subsampling problem. 

iMarconi et alJ l2003at) have shown how to derive the expressions used 
to compute the model quantities which must be compared with the 
observed ones. We summarize here their results and present a new 
way to compute model kinematical quantities which overcomes the 
"subsampling" problem described below. 

Consider a reference frame xy on the plane of the sky with the 
slit aligned along the y axis. The slit center has coordinate x = x and 



the points within the slit are characterized by Xq — Ax < x < Xq + Ax 
where 2Ax is the slit width projected onto the plane of the sky. The 
spectrograph provides matrixes 5 y where ( represents the wavelength 
direction and j represent the direction along the slit. For ( = 1, . . . , N 
(where N is the number of pixels along dispersion), the row Sy = 
*P/(Wj) is the spectrum extracted at a given position y>j along the slit 
and Wj is the velocity at pixel i. Such spectrum is extracted from an 
aperture which, on the plane of the sky, has coordinates x - Ax < x < 
x - Ax and y, - Ay < y < yj — Ay, where 2Ay is the pixel size along 
the slit. 

The kinematical quantities can be computed as: 



IMarconi etafl <2003ah show that: 



Vj) = 


J y ¥ j (w)dw = 2 Aw J dxj dy JJ dx'dy' P 


(Vlj) = 


^ w fj(yv) dw = 








= 2Aw 1 dx 1 dy II dx'dy' 


WoP 




(v 2 Ij) = 


J w 2 <P/w) dw = 








= 2 Aw f dx f dy ff dx'dy' 

sJx\ Jyi -J-J -co 


, Aw , 
n + — +o- 


9 



where x\?_ = xq ± Ax, yi^ = yi ± Ay, wo = v(x',y') + k(x - xq) and 
f = I(x',y')P(x - x', y - y'). v(x',y') is the velocity along the line of 
sight due to disk rotation, fc (x-Xp) is th e spurio us velocit y due to light 
entering off the slit center (Macieiewski & Binnev 2001) and P is the 
Point Spread Function (PSF). 2Aw is the width of the pixel in velocity 
andcri s the intrinsic velocity dispersion of the disk. See M arconi et alJ 
1 2003a) for more details and the derivation of these formulas. 

In order to compute (/,), {vlj} and (v 2 Ij), the more direct way is 
to build regular grid in the x',y' plane and do the convolution with 
the PSF b y using, for instance, the Fast Fourier Transform algorith 
(e.g. iPress et"ai]|l992h . After convolution, (/,•), {vlj} and (v 2 Ij) can 
then be computed by integrating the grids over the apertures from 
which observed spectra have been extracted. The requirement of reg- 
ular grids is imposed by the FFT algorithm and the grid should be fine 
enough to accurately sample the PSF; for instance, with a gaussian 
PSF, grids should have sampling steps less than or equal to 1/3 of the 
gaussian tr. The grid must also be large enough that edge effects are 
negligible over the aperture areas where the final integration is per- 
formed; for instance, with a gaussian PSF, it is enough to extend the 
grid by 5cr on all sides. 

There is a problem of subsampling when the intrinsic surface 
brightness distribution / is characterized by spatial scales which are 
much smaller that those of the PSF. A fine sampling of / results then 
in grids which are very large (more than, eg, 1024 x 1024 points) and 
which make computational times too long. However, while it is very 
important to finely sample / where there are the largest velocity gra- 
dients close to the BH location, this is a waste of computational re- 
sources where the velocity is slowly varying. A possibility to solve 
this problem is to invert the order of integration in the above expres- 
sions, i.e. first compute the integration over the apertures and then 
compute the FFT. In this way the small scales in / are removed by the 
integration process and the sampling must be fine only with respect to 
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Fig. A.l. Top: instrument convolved kinematical quantities computed with different grid sampling. "New" indicates the compu- 
tations with the new method described in the text. Bottom: as in the top panel but with the BH position offset by 0.02", 0.02" 
with respect to the peak of the surface brightness distribution. 



the PSF. The formulas to be used after the inversion of the integrals 
are trivially computed. Defining 

7(x,y) = I(x,y) 

v(x) = k * (x - x Q ) 

vl(x,y) = v(x,y)/(x,y) +7(x,y)v(x) 

v 2 I(x, y) = v 2 (x, y)I{x, y) + 7(x, y)v(x) 2 + 2 * vl(x, y)v(x) 



dx dy/(x,y) 

Jxq-Ax Jyi—Ay 
pxy+Ax ryj+Ay 

I dx I dy vl(x,y) 

Jxn-Ax Jy,- 



and 

I(.Xo,yj) 
vl(x ,yj) 

r-xo+Ax nyj+Ay 

v 2 /(x ,y ; ) = I dx I dyv 2 /(x,y) 

sJ x^-Ax Jyj-Ay 

the final expressions can be written as: 

*¥j(yv) dw = 2Aw I I dx'dy"Pl(x ,yj) 

CO \J\J — CO 



Jx -Ax Jyj-Ay 
r*xo +Ax r*y j + Ay 



<v/ 



X+co 
CO 



= 2AwJjf dx'dy'f [v/(x , y,) - v(x) * /(x„, y,)] 
<v 2 /;) = f w 2, F J -(w)dw = 2Awff dx'dy'!Ph77(xo,y ; ) 

*J —CO \_J \_) — CO 

-2 * v(x) * v/(x ,y,) + v(x) 2 * ZOo.y,)] 



where !P = P(jco — x 1 ,yj ->'')• The integrals over the apertures can be 
computed, e.g., with adaptive grids in order to reduce computational 
times. 

We now show the effects of subsampling on model values and 
how the above way of computation can solve them. We consider an 
intrinsic surface brightness distribution which is circularly symmetric 
on the plane of the sky and is described by a sum of two exponentials 

as 



I(r) = /„ eT rlra + Zj e^' 



(A.l) 



(w) dw : 



where r is the distance from the center and we assume the follow- 
ing parameter values: IJI = 0.01, r = 0.05", and r x = 1.00". The 
'observed' system has the same parameters as Centaurus A except 
for M BH = 10 8 M Q , M/L=l, and the gas disk is inclined by 45 deg with 
respect to the line of sight. The rotating gas disk, which has an intrin- 
sic constant velocity dispersion of 70 km s -1 , is observed with ISAAC 
with 0.5" seeing (FWHM). The slit is placed along the major axis of 
the disk. 
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The effects of subsampling are shown in Fig. IA, lh . The solid 
line represent the model values computed using the method described 
above, while the other lines represents computations with the com- 
monly adopted method and with different sampling steps. It is clear 
that the observed surface brightness distribution and velocity are very 
little affected by the subsampling problem. However the velocity dis- 
persion can be underestimated even by a factor 1.5 for the largest 
sampling used which is nonetheless a factor 5 smaller than the PSF 
o\ The effects on velocity and velocity dispersion can also be seen 
if the BH location and peak of the surface brightness distribution are 
not located on one of the grid points. In Fig. lA. lb the BH is offset by 
only 0.02", 0.02" along and across the slit with respect to the center 
of surface brightness distribution. In this case all gas velocity disper- 
sions are close to the correct value indicated by the solid line but the 
velocity values are starting to be affected by the sampling problem. 



